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Two  synthetic  light  curves  computed  from  theoretical 
astrophysical  models  of  W  Ursae  Majoris  systems  are  discussed. 
Solutions  of  these  light  curves,  based  on  the  geometrical 
Russell  model  and  the  Russell-Merrill  method  of  solution  of 
the  light  curves  of  eclipsing  binaries,  are  presented.   The 
relatively  shallow  minima  caused  by  the  partial  eclipses 
of  the  first  synthetic  light  curve  lead  to  a  problem  of 
indeterminacy  in  the  solution.   The  "observed"  points  in  the 
shoulders  of  both  synthetic  light  curves  fall  below  the 
theoretical  light  curves  predicted  by  the  solutions.   The 
addition  of  third  light  in  the  solution  of  the  two  synthetic 


light  curves  improves  the  fit  of  the  solutions  to  the  light 
curves;  however,  there  is  no  sound  basis  for  adding  this 
third  light.   The  orbital  elements  predicted  by  the  Russell- 
Merrill  solutions  of  the  two  synthetic  light  curves  are  not 
at  all  close  to  the  orbital  elements  used  to  generate  these 
light  curves  from  the  theoretical  astrophysical  models.   In 
particular  the  Russell-Merrill  solution  underestimates  the 
sizes  of  the  components.   It  is  concluded  that  the  Russell 
model  is  not  compatible  with  the  theoretical  astrophysical 
models  used  to  generate  the  synthetic  light  curves. 

Observations  of  the  system  of  44i  Bootis ,  an  example 
of  a  W  Ursae  Majoris  system,  are  also  discussed.   The  times 
of  minimum  light  indicate  that  an  increase  in  the  period 
occurred  in  1967. 


CHAPTER  I 


INTRODUCTION 


The  W  Ursae  Majoris  systems  are  eclipsing  variable 
stars  whose  light  curves  have  maxima  which  are  strongly 
curved  and  minima  which  are  nearly  equal  in  depth.   These 
systems  have  periods  which  average  approximately  one-half 
day.   Orbital  solutions  in  the  literature  indicate  that 
W  Ursae  Majoris  systems  are  close  binaries  whose  separations 
are  less  than  the  dimensions  of  the  components.   Many 
complexities  which  are  caused  by  the  proximity  of  the  com- 
ponents appear  in  the  light  curves  (1) .   Most  of  the  solu- 
tions of  the  light  curves  of  W  Ursae  Majoris  systems  in 
the  literature  are  based  upon  a  geometrical  model  first 
proposed  by  Russell  (2) .   Lucy  has  proposed  that  some  or 
most  of  the  W  Ursae  Majoris  systems  may  be  true  contact  sys- 
tems, whose  common  boundary  follows  a  single  equipotential 
surface  (3).   It  is  the  purpose  of  this  thesis  to  apply  the 
Russell-Merrill  method  of  solution  of  the  light  curves  of 
eclipsing  binaries  (2,4,5)  to  synthetic  light  curves  from 
astrophysical  models  similar  to  Lucy's  model  in  order  to 
determine  whether  the  geometrical  Russell  model  is  compati- 


ble  with  the  astrophysical  models.   Lucy  (3),  Rucinski  (6), 
and  Mochnacki  and  Doughty  (7)  have  written  computer  programs 
for  computing  theoretical  light  curves  of  W  Ursae  Majoris 
system;  Mochnacki  and  Doughty  have  published  a  trial  and 
error  fit  to  the  system  AW  Ursae  Majoris  using  their  program. 
Wilson  and  Devinney  (8)  have  published  a  general  procedure 
for  computing  light  curves  of  close  binaries  which  includes 
the  W  Ursae  Majoris  systems  as  a  special  case.   This  pro- 
cedure is  now  being  applied  to  selected  contact  systems  and 
the  results  will  be  published  soon.   Two  synthetic  light 
curves  have  been  computed  for  this  study,  one  by  Rucinski  (6) 
and  one  by  Wilson  (9) .   In  addition,  observations  of  the 
system  of  44i  Bootis,  an  example  of  a  W  Ursae  Majoris  system, 
are  discussed. 

The  Russell  Model 

The  Russell  model  and  the  Russell-Merrill  method  of 
solution  of  the  light  curves  of  eclipsing  binaries  are 
discussed  in  detail  by  Russell  and  Merrill  (2) .   A  brief 
discussion  is  given  here  in  order  to  define  the  notation 
used.   The  Merrill  tables  and  nomographs  for  solution  of 
light  curves  of  eclipsing  binaries  are  based  upon  a  spherical 
model  (4,5).   It  is  assumed  that  the  orbit  is  circular, 
that  the  stars  are  spherical,  and  that  they  appear  darkened 
at  the  limb  by  a  linear  cosine  darkening  law.   This  limb 


darkening  (x)  may  differ  for  the  two  components.   The  light 
of  each  star  is  constant  for  the  spherical  model;  therefore 
the  light  of  the  system  outside  of  eclipse  is  also  constant. 
It  is  convenient  to  take  the  light  of  the  system  outside 
of  eclipse  as  the  unit  of  light  and  the  radius  of  the  spheri- 
cal orbit  as  the  unit  lof  component  dimensions.   Then  the 
components  are  defined  by  the  following  parameters: 

Radius  of  the  larger  (greater)  star   r 

y 

Radius  of  the  smaller  star  r 

Inclination  of  the  orbit  i 

Light  of  the  larger  (greater)  star     L 

y 

Light  of  the  smaller  star  L 

with  L  +  L  =1.   If  9  is  the  longitude  in  orbit  (from 
g    s 

conjunction) ,  then  the  apparent  distance  between  the  centers 
of  the  disks  is  given  by 

2      2       2      2 

6   =  cos  i  +  sin  i  sin  6 


Setting  p  =  (6  -  r  )/r  ,  the  eclipse  will  be  absent,  partial, 
g    s 

or  complete,  for  p>l,  l>p>-l,  orp<  -1  respectively. 

Setting  k  =  r  /  r   then  6  =  r  (1  +  kp) .   The  quantities  k 
^      s   g  g 

and  p  are  dimensionless,  and  their  values  completely  define 

the  geometrical  circumstances  of  a  given  phase. 

If  f   and  f  represent  the  fractions  of  the  light  of 
g      s 

the  two  stars  which  are  obscured  at  any  phase  of  the  eclipse 
of  either,  and  i   is  the  normalized  value  of  the  light  re- 


ceived  from  the  whole  system: 


«  =  L^d  -  f^)  +  L^d  -  f^)  =  1  -  Lgf^-  L^f^ 


For  tabular  purposes  Russell  and  Merrill  express  these  in 
terms  of  two  other  functions  a  and  x,  where  a  is  the  ratio 
of  light  lost  at  any  phase  of  an  eclipse  to  the  loss  at 
internal  tangency,  and  x  is  that  of  the  latter  to  the  whole 
light  of  the  star.   Then  for  the  light  at  any  phase  during 
an  occultation  (the  larger  star  in  front) 


1  -  £   =  f  Lg  =  Lgtt   (Xg,k,p) 


and  during  a  transit  (the  smaller  star  in  front) 


1  _  5,tr  ^  ftr^^  ^  Lg  x(Xg,k)a  ^^(Xg,k,p) 


It  is  convenient  to  use  the  x-functions  of  Russell 
and  Merrill  to  determine  the  solution  of  the  light  curve. 
Defining  n  =  a/a^  where  the  zero  subscript  refers  to  the 
value  of  the  parameter  at  mid-eclipse,  the  X'f unctions  are 
given  as 


X(x,k,a^;n)=   ^^^^^in)        ^      [l+kp(x,k,nao)]^ 
sin2e(n=0.5)    [1+kp (x,k,0 .5a^) ] ^ 


_  [l+kp(x,k/ao)]^ 
[l+kp(x,k,a^)]2 


These  x~functions  have  been  tabulated  by  Merrill  (4). 

Russell  and  Merrill  have  shown  that  a  system  consisting 

of  two  similar  triaxial  ellipsoids  with  semi-axes  a  ,  b  , 

g   g 

c  and  a  =  ka  ,  b  =  kb  ,  c  =  kc   can  be  rectified  to 
g      s     g'   s     g    s     g 

the  spherical  model  making,  certain  approximations  which 
involve  the  gravity  and  reflection  effects.   The  a-axis  of 
each  ellipsoid  is  along  the  line  joining  the  centers  of  the 
components,  the  c-axis  is  parallel  to  the  axis  rotation  of 
the  system,  and  the  b-axis  is  in  the  third  mutually  perpen- 
dicular direction.   A  mean  radius  r  =  (a  +  b  +  c)/3  may 
also  be  defined  for  these  ellipsoids.   The  fundamental  geo- 
metrical equations  for  this  model  (i.e.   the  Russell  model) 
may  be  written  as 


2        2      2       2        2 

cos  i  +  sin  x  sin  9  =  a   (1  +  kp) 
r        r         g 


where  i   and  0  are  the  rectified  inclination  and  phase  angle 

(orbital  longitude)  respectively.   This  equation  is  identical 

in  form  with  the  equation  for  spherical  stars.   Thus  the 

observed  intensity  and  phase  angle  can  be  rectified  in  such 

a  way  as  to  produce  a  rectified  light  curve  which  will  be 

nearly  that  produced  by  the  eclipse  of  a  pair  of  spherical 

stars  of  radii  a  ,  a  =  ka  ,  and  inclination  i„. 
g   s     g  r 


CHAPTER  II 
FIRST  SYNTHETIC  LIGHT  CURVE 

A  synthetic  light  curve  for  a  W  Ursae  Majoris  type 
eclipsing  binary  was  generated  by  Dr.  S.  M.  Rucinski  from 
Lucy's  model  (6)  .  The  "observational"  data  for  this  light 
curve  consisted  of  37  values  of  the  normalized  light  (or 
flux)  as  a  function  of  phase  angle,  with  the  phases  given 
in  5°    intervals  from  0°  to  180°.   The  other  half  of  the 
light  curve  (180°  to  360°)  was  assvimed  to  be  symmetrical. 
An  additional  ten  points  were  generated  later  to  define 
better  the  centers  of  the  eclipse  regions  of  the  light 
curve.   This  "observational"  data  is  listed  in  Table  I 
and  plotted  in  Figure  I.   This  synthetic  light  curve  was 
to  be  treated  as  observational  data  and  solved  by  the 
standard  Russell-Merrill  method  of  solution  of  light  curves 
of  eclipsing  binaries .   No  additional  information  about  the 
nature  of  the  derivation  of  this  light  curve  was  to  be  used 
in  the  solution. 

Rectification 
The  first  step  in  the  process  of  getting  a  solution 


TABLE  I 


Observations  of  First  Synthetic  Light  Curve 


PHASE 


INTENSITY 


PHASE 


INTENSITY 


0.00 

0.45242 

95.00 

0.99343 

2.50 

0.44826 

100.00 

0.98470 

5.00 

0.45216 

105.00 

0.96954 

7.50 

0.46429 

110.00 

0.94766 

10.00 

0.48257 

115.00 

0.92548 

12.50 

0.49728 

120.00 

0.90144 

15.00 

0.51875 

125.00 

0.86813 

17.50 

0.54518 

130.00 

0.83688 

20.00 

0.56516 

135.00 

0.80528 

22.50 

0.58812 

140.00 

0.77124 

25.00 

0.61087 

145.00 

0.73817 

30.00 

0.65572 

150.00 

0.70216 

35.00 

0.70708 

155.00 

0.66023 

40.00 

0.74969 

157.50 

0.64105 

45.00 

0.79267 

160.00 

0.62372 

50.00 

0.83087 

162.50 

0.60589 

55.00 

0.86916 

165.00 

0.58505 

60.00 

0.90446 

167.50 

0.56515 

65.00 

0.92940 

170.00 

0.55371 

70.00 

0.95640 

172.50 

0.54237 

75.00 

0.97754 

175.00 

0.53530 

80.00 

0.99117 

177.50 

0.53398 

85.00 

0.99632 

180.00 

0.53866 

90.00 

0.10000 

for  the  light  curve  was  an  analysis  of  the  light  outside 

the  eclipses  in  order  to  arrive  at  a  rectification  of  the 

light  curve  to  the  spherical  model.   Two  methods  were  used 

for  this  analysis:   Merrill's  graphical  method  (10)  and 

a  least  squares  Fourier  analysis  of  the  material  outside  of 

eclipse. 

Let  the  light  outside  eclipse  be  represented  by  a 
truncated  Fourier  series  of  the  form: 


I  =  A   +  A,  cos  0  +  A^  cos  29  +  A-,  cos  36  +  A.  cos  4( 
o     1  2  3  4 


(Since  the  light  curve  is  symmetrical  about  9  =  180°,  it 
is  not  necessary  to  include  sine  terms  in  the  above  Fourier 
series  i')   Following  Merrill's  graphical  method,  a  and  b 
represent  readings  for  9  and  180"  -  6  on  the  light  curve. 
(The  given  "observational"  points  were  used  rather  than 
reading  from  a   freehand  curve  since  the  scatter  of  the 
points  was  small.)   It  immediately  follows  that 


2-(a  +  b)  =  A^  +  A^   cos  29  +  A^  cos  4( 


p-(a  -  b)  =  A-  cos9  +  A^  cos  3( 


Letting  C,  =  ^(a  -  b)  and  C^  =  yC^  +  b) ,  by  simple  trigo- 
nometric substitution, 

C,  =  (A,  -  3A^)    cos  0  +  4A2  cos  36 


10 


C2  =  (Aq  -  A^)  +  A2  cos  26  +  2A^  cos^  29 


Therefore  a  plot  of  C^^  versus  cos 6  would  have  the  form  of 
a  cubic  and  a  plot  of  Q.^   versus  cos  20  would  have  the  form 
of  a  parabola.   Such  plots  are  given  in  Figure  2.   The  plot 
of  C^  versus  cos  9  shows  the  presence  of  a  considerable  A- 
(cos  30)  term  and  indicates  the  possible  presence  of  higher 
order  odd  cosine  terms.   The  plot  of  Q^   versus  cos  20  is 
essentially  linear,  indicating  that  the  A^   term  is  negli- 
gible.  The  plotted  points  seem  to  fall  below  this  straight 
line  somewhere  around  40°,  thus  locating  0   (the  value  of 
9  at  external  tangency)  to  a  first  approximation.   Values 
of  the  Fourier  cosine  coefficients  were  then  derived  from 
the  plots  and  are  listed  in  Table  II. 

Because  of  the  similarity  of  the  "colors"  of  the  two 
"stars"  and  the  small  difference  in  the  depths  of  the  two 
minima,  rectification  coefficients  for  the  reflection  effect 
were  obtained  in  the  following  manner  (2): 


^1  =  -^1 


C   =  0.090  sin^  9 
o  e 


C,  =  0.030  sin^  0 
2  e 


The  eclipses  were  assumed  to  be  partial  and  values  were 
derived  for  0^  =  39°;  these  values  are  listed  in  Table  II, 
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TABLE  II 

Graphical  Rectification  Coefficients  for 
First  Synthetic  Light  Curve 


e 

a 

b 

^1 

^2 

0° 

.4534 

.5387 

-.0432 

.4956 

5 

.4522 

.5353 

-.0416 

.4938 

10 

.4826 

.5357 

-.0356 

.5182 

15 

.5188 

.5851 

-.0332 

.5520 

20 

.5652 

.6237 

-.0293 

.5945 

25 

.6109 

.6602 

-.0247 

.6356 

30 

.6557 

.7022 

-.0233 

.6790 

35 

.7051 

.7382 

-.0166 

.7217 

40 

.7497 

.7712 

-.0108 

.7605 

45 

.7927 

.8053 

-.0063 

.7990 

50 

.8309 

.8369 

-.0030 

.8339 

55 

.8692 

.8681 

.0006 

.8687 

60 

.9045 

.9014 

.0016 

.9030 

65 

.9294 

.9355 

.0020 

.9275 

70 

.9564 

.9477 

.0044 

.9521 

75 

.9775 

.9695 

.0040 

.9735 

80 

.9912 

.9847 

.0033 

.9880 

85 

.9963 

.9934 

.0015 

.9949 

90 

1 

1 

0 

1 

^o 

=   0.7990 

^1 

=  -0.0232 

^2 

=  -0.2010 

^3 

=  -0.0147 

^4 

=   0 

C   =  0.0354 
o 

^1  =  -^1 
C^   =  0.0118 


0.3353  0.6 

z  =   0.2724   for  x  =   0.8 

0.2180  1.0 
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Also  listed  in  Table  II  is  the  value  of  z,  given  by 


-4(A2  -  C2) 


N(^o  -  ^o  -  ^2  -^  ^2^ 


where  N  depends  on  the  assiomed  limb  darkening  (x)  ',   values 
were  adopted  such  that  N  =  2.6,   3.2,  or  4.0  when  x  is 
assxuned  to  be  0.6,  0.8  or  1.0  respectively. 

Fourier  coefficients  were  also  computed  by  a  least 
squares  Fourier  analysis  of  the  light  outside  of  eclipses 
(45°  to  135°) .   The  results  of  this  analysis  are  listed  in 
Table  III. 

Both  sets  of  rectification  coefficients  were  then 
used  to  compute  a  rectified  light  curve.   The  rectified 
intensity  is  given  by 


I+C  +C-,    cos  e+C^  cos  29 -A^  cos  36 -A^  cos  4( 
I"  = 


^^  ^    ,^.^^  ^^^    .  .V.2  —  ^.  «3  v,^^  ...  «4 


A  +C  +  (A^+Co)  cos  2( 
00     2   2 


and  the  rectified  phase  by 


sin^e  =  -^ 


■I       2. 
1-z  cos  f 


The  rectified  points  are  listed  in  Tables  IV  and  V 
and  plotted  in  Figures  3  and  4. 

Several  things  are  apparent  from  an  examination  of 
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TABLE  III 


Least  Squares  Rectification  Coefficients  for 
First  Synthetic  Light  Curve 


Aq 

=   0.79532 

+ 

0.00123 

^1 

=  -0.02208 

+ 

0.00114 

^2 

=  -0.20699 

+ 

0.00195 

^3 

=  -0.01304 

+ 

0.00066 

^4 

=  -0.00350 

+ 

0.00090 

OBSERVED 

THEORETICAL 

0-C 

e 

INTENSITY 

INTENSITY 

0.7927 

0.7924 

0.0003 

45.0000 

0.8309 

0.8317 

-0.0008 

50.0000 

0.8692 

0.8687 

0.0005 

55.0000 

0.9045 

0.9026 

0.0019 

60.0000 

0.9294 

0.9322 

-0.0028 

65.0000 

0.9564 

0.9570 

-0.0006 

70.0000 

0.9775 

0.9763 

0.0012 

75.0000 

0.9912 

0.9898 

0.0014 

80.0000 

0.9963 

0.9973 

-0.0010 

85.0000 

1.0000 

0.9988 

0.0012 

90.0000 

0.9934 

0.9944 

-0.0010 

95.0000 

0.9847 

0.9845 

0.0002 

100.0000 

0.9695 

0.9693 

0.0002 

105i0000 

0.9477 

0.9495 

-0.0018 

110.0000 

0.9255 

0.9257 

-0.0002 

115.0000 

0.9014 

0.8986 

0.0028 

120.0000 

0.8681 

0.8689 

-0.0008 

125.0000 

0.8369 

0.8375 

-0.0006 

130.0000 

0.8053 

0.8052 

0.0001 

135.0000 
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TABLE  IV 

Rectified  First  Synthetic  Light  Curve 
Using  Graphical  Rectification  Coefficients 


0(x=.6) 


0(x=.8) 


0(x=l) 


I"  (  e ) 


I"  (180+9 P; 


0.000 

3.066 

2. 

931 

6.125 

5 

857 

9.173 

8 

775 

12.204 

11 

680 

15.212 

14 

570 

18.194 

17 

439 

21.143 

20 

286 

24.058 

23 

108 

26.933 

25 

901 

29.768 

28 

665 

35.405 

34 

092 

40.658 

39 

.382 

45.825 

44 

.540 

50.810 

49 

.536 

55.573 

54 

.406 

60.279 

59 

.152 

64.793 

63 

.781 

65.000 

70.000 

75.000 

80.000 

85.000 

90.000 

2.826 
5.651 
8.468 
11.277 
14.074 
16.857 
19.624 
22.371 
25.099 
27.803 
33.141 
38.373 
43.498 
48.513 
53.423 
58.234 
62.953 


0.8331 
0.8254 
0.8278 
0.8404 
0.8599 
0.8714 
0.8906 
0 . 9145 
0.9264 
0.9405 
0.9522 
0.9693 
0.9870 
0.9950 
0.9996 
0.9990 
1.0006 
1.0020 
0.9967 
0.9986 
1.0001 
1.0000 
0.9976 
1.0000 


0.8493 
0.8413 
0.8410 
0.8480 
0.8599 
0.8703 
0.8916 
0.9124 
0.9268 
0.9388 
0.9517 
0.9778 
0.9905 
0.9959 
1.0003 
1.0009 
1.0014 
1.0054 
1.0018 
0.9995 
1.0009 
1.0001 
0.9983 
1.0000 
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TABLE  V 


Rectified  First  Synthetic  Light  Curve  Using  Least 
Squares  Rectification  Coefficients 


0.4524 
0.4483 
0.4522 
0.4643 
0.4826 
0.4973 
0.5188 
0.5452 
0.5652 
0.5881 
0.6109 
0.6557 
0.7051 
0.7497 
0.7927 
0.8309 
0.8692 
0.9045 
0.9294 
0.9564 
0.9775 
0.9912 
0.9963 
1.0000 
0.9934 
0.9847 
0.9695 
0.9477 
0.9255 
0.9014 
0.8681 
0.8369 
0.8053 


I" 


0.8469 

0.000 

0.8391 

2.945 

0.8412 

5.885 

0.8535 

8.816 

0.8726 

11.735 

0.8834 

14.637 

0.9019 

17.519 

0.9250 

20.377 

0.9359 

23.208 

0.9488 

26.011 

0.9595 

28.781 

0.9745 

34.221 

0.9904 

39.518 

0.9968 

44.668 

1.0004 

49.673 

0.9992 

54.538 

1.0006 

59.273 

1.0021 

63.891 

0.9971 

68.405 

0.9994 

72.831 

1.0012 

77.186 

1.0014 

81.488 

0.9990 

85.753 

1.0012 

90.000 

0.9990 

94.248 

1.0003 

98.513 

1.0002 

102.815 

0.9982 

107.170 

0.9998 

111.596 

1.0031 

116.110 

0.9992 

120.727 

0.9994 

125.463 

1.0002 

130.328 
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Table  V  continued. 


0.7712 
0.7382 
0.7022 
0.6602 
0.6411 
0.6237 
0.6059 
0.5851 
0.5652 
0.5537 
0.5424 
0.5353 
0.5340 
0.5387 


I" 


0.9977 

135.332 

0.9952 

140.483 

0.9857 

145.780 

0.9631 

151.220 

0.9521 

153.990 

0.9417 

156.792 

0.9290 

159.624 

0.9097 

162.482 

0.8897 

165.364 

0.8804 

168.266 

0.8694 

171.185 

0.8630 

174.116 

0.8639 

177.056 

0.8722 

180.000 
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either  of  the  rectified  light  curves.   The  first  is  that 
the  amount  of  scatter  of  the  data  points  in  eclipse  is  much 
larger  than  was  originally  expected.   This  scatter  will  be 
treated  as  "observational"  error  for  the  present  and  no  at- 
tempt will  be  made  to  explain  it  either  in  terms  of  the 
model  or  the  method  of  computation  of  the  data.   Another 
obvious  feature  of  the  rectified  minima  is  the  "brightening" 
at  the  centers  of  the  eclipses.   Again  no  attempt  will  be 
made  at  present  to  explain  this  effect;  however,  its  pres- 
ence^ creates  a  serious  uncertainty  in  both  the  depths  of 
the  eclipses  and  the  credibility  of  the  points  near  the 
center  in  terms  of  the  Russell  model. 

Solutions  from  Graphical  Rectification 

The  rectified  light  curve  produced  by  applying  the 
rectification  coefficients  derived  by  the  graphical  method 
was  used  for  the  first  attempts  at  a  solution.   The  phase 
was  rectified  with  a  z  based  on  N  =  2.6  corresponding  to  a 
darkening  x  =  0.6,  and  the  first  solutions  were  tried  on  the 
X  =  0.6  nomograph.   The  depths  of  the  eclipses  were  chosen 
as  1  -  ij^^   =0.18  and  1  -  l^^^^   =0.16  and  the  following 
values  were  read  off  the  plot  of  the  rectified  light  curve: 

n      0       sin^e     x^""       e     sin^e     x^®"" 

0.2     33?0       .29663      2.268      32?6     .29027    1.901 
0.5     21.2       .13077       1       23.0     .15267     1 
0.8     12.2       .044658      .342      15.0     .066987    .439 
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By  choosing  values  of   0(n  =  .5)  and   9{n  =  .8) 

PIT 

at  9  +  0?5,  a  permissible  range  of  .301  <_  %      (n  =  .8)    <_ 
.387  and  .394  <  ^^^^(n  =  .8)  <  .488  was  found.   Since 
^secjj^  =  .8)  >  x^^^^  ~  '^^  '    ^^   appears  that  the  primary 
eclipse  is  a  transit  and  the  secondary  eclipse  is  an  occul- 
tation.   Then  the  values  needed  for  the  depth  line  on  the 
nomograph  are 

(1  -  Jl°°)  +  (1  -  ll"^)    =0.34 

1  -  il^/i°^   =  0.2142 

1  -  £oVil°''  +  (1/50)  1  -  £°''/ii°°=  0.2180 


Taking  these  values  and  the  x(n  =  0.8)  values  given  above 
to  the  X  =  0.6  nomograph,  an  intersection  of  the  depth  line 
and  the  permissible  values  of  both  x  contours  was  found 
with  k  =  0.45,  p  =-.98,  a^^  =  .9987,  a^^  =  .9928.   When 
these  elements  are  taken  to  the  x  tables,  they  produce  the 
following  points  on  the  light  curve: 


X°^ 


0.0  2.823 

0.2  1.725 

0.5  1 

0.8  .448 


Comparing  these  values  with  those  taken  from  the 
"observed"  light  curve,  several  things  are  apparent.   First, 


•  2„ 
sin  0 

0 

x"-^ 

2 

sxn  0 

0 

.44122 

41?6 

3.374 

.44122 

41°6 

.26961 

31.3 

1.913 

.25016 

30.0 

.15629 

23.3 

1 

.13077 

21.2 

.070020 

15.3 

.376 

.049170 

12.8 

24 


the  fit  from  the  half-way  point  down  is  moderately  satis- 
factory for  this  preliminary  stage  of  solution.   Second, 
the  fit  at  the  shoulders  is  very  bad;  the  computed  curve 
is  much  narrower  than  the  "observed"  curve.   The  x (n  <  .5) 
values  for  both  eclipses  need  to  be  increased  a  significant 
amount  while  making  only  small  changes  in  the  x(n  >  .5) 
values,  in  order  to  fit  the  "observed" light  curve  with  the 
chosen  depths  of  the  eclipses. 

Since  the  nomographic  solution  on  darkening 
X  =  0.6  as  described  above  was  not  satisfactory,  other 
possibilities  were  explored.   The  first  approach  was  to 
try  nomographic  solutions  with  other  values  of  darkening. 
Table  VI  summarizes  the  results  of  these  attempts.   The 
main  conclusion  from  this  exploration  of  solutions  with 
different  values  of  darkening  is  that  the  fit  will  improve 
with  increasing  darkening.   There  are  two  reasons  for  this 
improvement.   First  the  x (n  <  .5)  values  for  the  transit 
eclipse  tend  to  increase  with  increasing  darkening,  and 
second,  the  value  of  z  is  smaller  for  larger  darkening. 
With  a  smaller  z,  the  rectification  of  the  phase  tends  to 
make  the  shoulders  narrower  relative  to  the  half-width 
than  with  a  larger  z.   Thus,  it  appeared  that  a  darkening 
of  X  =  .8  or  X  =  1  should  be  used  for  further  trial 
solutions.   The  results  of  these  solutions  with  k  =  .45 
(see  Tables  VII  and  VIII  and  Figures  5  and  6)  were  not  at 
all  satisfactory.   It  was  still  not  possible  to  fit  the 
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TABLE  VI 
Solutions  with  Different  Values  of  Liitib  Darkening 


X 

0.2 

0.4 

0.6 

0.8 

1.0 

k 

0.4625 

0.4625 

0.45 

0.45 

0.45 

Po 

-0.93 

-0.90 

-0.98 

-0.94 

-0.90 

a^oc 

0.9868 

0.9799 

0.9987 

0.9943 

0.9905 

a  tr 
o 

0.9809 

0.9627 

0.9928 

0.9698 

0.9375 

X^^{n= 

.8) 

0.401 

0.405 

0.448 

0.452 

0.452 

X^^(n= 

.8) 

0.383 

0.370 

0.376 

0.361 

0.351 

X°^(n= 

.2) 

1.812 

1.799 

1.725 

1.712 

1.703 

X^^(n= 

.2) 

1.867 

1.908 

1.913 

1.957 

1.991 
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TABLE  VII 

Solution  for  x  =  0.8,  1  -  il^^  =  0.18,  1  -  Jl^®^=  0.16, 
k  =  0.45,  p  =  -0.94,  a°^  =  0.9943,  a^^  =  0.9698 


n 

X°" 

2 

sin  9 

0 

tr 

X 

sin^G 

0 

0.0 

2.827 

.47732 

4327 

3.650 

.47732 

43?7 

0.1 

2.047 

.34562 

36.0 

2.464 

.32222 

34.6 

0.2 

1.712 

.28905 

32.5 

1.957 

.25592 

30.4 

0.4 

1.203 

.20312 

26.8 

1.265 

.16542 

24.0 

0.5 

1 

.16884 

24.3 

1 

.13077 

21.2 

0.6 

.817 

.13794 

21.8 

;766 

.100170 

18.5 

0.8 

.452 

.076316 

16.0 

.361 

.047208 

12.6 

0.9 

.276 

.044600 

12.2 

.181 

.023669 

8.9 

0.95 

.1734 

.029277 

9.9 

.0925 

.012096 

6.3 
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TABLE  VIII 


Solution  for  x  =  1.0,  1  -  l^^  =   0.18,  1  -  £^®^  =0.16 

'      o  o 

k  =  0.45,  p  =  -0.90,  a°^  =  0.9905,  a^^  =  0.9375 
'^O  o  o 


0.0  3.801   2.911  .47732   43?7  5.412   3.939  .47732   43.7 

0.1  2.257   2.054  .33675   35.5  3.082   2.532  .30687   33.6 

0.2  1.625   1.703  .27922   31.9  2.186   1.991  .24132   29.4 

0.4     .7274  1.205  .19750   26.4    .9965  1.273  .15430   23.1 

0.5     .3590    1  .16396   23.9    .5439    1  .12119   20.4 

0.6     .0191   .811  .13302   21.4    .1465   .760  .092118  17.6 

0.8  -  .6275   .452  .074151  15.8  -.5306   .351  .042584  11.9 

0.9  -  .9684   .263  .043116  12.0  -.8290   .171  .20755    8.3 

0.95  -1.1648   .1539  .025236   9.1  -.9717   .0851  .010315   5.8 
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shoulders  of  the  eclipses,  especially  the  primary  eclipse, 
while, at  the  same  time,  fitting  the  rest  of  the  eclipse 
curves.   In  particular,  with  k  =  .45,  the  solutions  for 
a  circular  orbit  were  too  wide  at  the  halfway  point  of 
the  secondary  eclipse  and  too  narrow  on  the  shoulders  of 
the  primary  eclipse,  with  lesser  problems  elsewhere. 

Because  of  the  problems  described  above,  it  was 
necessary  to  abandon  the  chosen  value  of  k  =  0.45  and  to 
explore  other  possibilities  of  k  on  the  x  =  0.6  nomograph, 
staying  on  the  depth  line  given  above.   Table  IX  summarizes 
the  results  of  these  explorations.   From  an  inspection  of 
this  table,  it  was  decided  to  attempt  trial  solutions  with 
k  in  the  range  .50  to  .60  for  darkenings  x  =  0.8  and  total 
darkening.   Some  of  these  solutions  are  given  in  Tables 
X,  XI,  and  XII.   While  these  solutions  are  an  improvement 
over  previous  ones,  the  major  problem  of  the  fit  of  the 
shoulders,  especially  in  the  primary  eclipse,  has  not  been 
alleviated. 

Because  of  the  recurrent  problem  with  the  fit  of 
the  shoulders  of  the  eclipses,  it  was  felt  that  a  different 
approach  might  prove  helpful.   This  approach  was  to  choose 
a  0 (n  =  .5)  for  the  primary  eclipse  and  a  9  .   From  these 

values,  X  (ii  -   0)'s  were  derived.   These  X's  were  taken  to 

oc 
the  X  =  0 . 8  tables  and  values  of  k  and  a   were  obtained 

o 

and  these  shape  curves  were  plotted  (see  Table  XIII  and 
Figure  7) .   The  depth  curve  derived  from 
a|°  =  (1  -  1°"")    +  (1  -  lT)/q^ 
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TABLE  IX 

Solutions  for  x  =  0.6,  1  -  l^^/n!^^   =0.2142 

o   o 

(1  -9.^)    +  (1  -  ii°^)  =0.34 


k  1.00     .85  .75     .60     .55     .50      .45 

Pq         .0956   .0000  -.1000  -.3538  -.4820  -.6789  -.9800 

a^oc       .3400   .4136  .4900   .6800   .7700   .8900    .9987 

a^tr       .3400   .3788  .4363   .6018   .6897   .8224    .9928 

X°^(n=.8)   .346    .346  .347    .354    .361    .379     .448 

x''^^(n=.8)   .346    .344  .343    .343    .344    .348     .376 

oc 

X   (n=.2)  1.969   1.963  1.956  1.927  1.902   1.853  1.725 

X^^(n=.2)  1.969      1.976  1.982  1.986  1.983      1.972  1.913 
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TABLE  X 

Solution  for  x  =  0.8,  l-;iP^  =  0.18,  1-  l^^^   =  0.16, 

o  o 

k  =  0.50,  p  =  -0.6608,  a°^  =  0.89,  a^^  =  0.7915 


n 

oc 

X 

sin  0 

G 

x^=^ 

2 

sin  0 

0;^ 

0.0 

3.217 

.47732 

43?7 

3.746 

.47732 

43°7 

0.1 

2.261 

.33547 

35.4 

2.524 

.32161 

34.6 

0.2 

1.838 

.27271 

31.5 

1.995 

.25421 

30.3 

0.4 

1.241 

.18413 

25.4 

1.277 

.16272 

23.8 

0.5 

1 

.14837 

22.7 

1 

.12742 

20.9 

0.6 

.781 

.115880 

19.9 

.758 

.096585 

18.1 

0.8 

.384 

.056976 

13.8 

.345 

.043960 

12.1 

0.9 

.195 

.028933 

9.8 

.166 

.021152 

8.4 

0.95 

.0992 

.014719 

7.0 

.0816 

.010398 

5.9 
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TABLE  XI 

Solution  for  x  =  0.8,  1  -  l^^  =   0.18,  1  -  l^^^=   0.16, 

o  o        ' 

k  =  0.55,  p^  =  -0.4798,  a°^  =  0.78,  a^^  =  0.6659 
'^o  '   o  o 


n-: 

X°^ 

sin^G 

0 

tr 

X 

sin^G 

G 

0.0 

3.382 

.47732 

43?7 

3.764 

.47732 

43?7 

0.1 

2.347 

.33124 

35.1 

2.535 

.32147 

34.5 

0.2 

1.892 

.26703 

31.1 

2.002 

.25388 

30.3 

0.4 

1.254 

.17698 

24.9 

1.278 

.15207 

23.7 

0.5 

1 

.141135 

22.1 

1 

.126812 

20.9 

0.6 

.771 

.108815 

19.3 

.756 

.095870 

18.0 

0.8 

.366 

.051655 

13.1 

.343 

.043497 

12.0 

0.9 

.180 

.025404 

9.2 

.164 

.020797 

8.3 

0.95 

.0893 

.012603 

6.4 

.0803 

.010183 

5.8 
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TABLE  XII 

Solution  for  X  =  1.0,  1  -  Ji^^  =  0.18,  1  -  H^^^   =  0.16, 

k  =  0.60,  p  =  -0.3225,  a°^  =  0.68,  a'^^  =  0.5332 
'  '^O               o  o 


, oc      oc  .  2n  r\           .  tr     tr  .  2^  „ 

n     t\)                 X  sin  0  0    ij;  X  sin  0  0 

0.0    4.495  3.582  .47732  43?7  6.349  3.982  .47732  43?7 

0.1    2.910  2.392  .31875  34.4  4.186  2.562  .30710  33.7 

0.2    2.276  1.916  .25532  30.3  3.349  2.012  .24117  29.4 

0.4    1.4016  1.259  .16777  24.2  2.233  1.280  .15343  23.1 

0.5    1.0562     1  .133255  21.4  1.8073    1  .119869  20.3 

0.6    .7471     .768  .102340  18.7  1.4350  .756  .090621  17.5 

0.8    .2032     .360  .047972  12.7    .0843  .342  .040995  11.7 

0.9   -.0424     .175  .023320  8.8    .5327  .163  .019539  8.0 

0.95  -;16b2     .0868  .011567  6.2    .4061  .0800  .009590  ::5;6 
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TABLE  XIII 


Shape  Curve  for  0  =  50°,  x°^  =  3.955, 


X^^  =  4.605,  X  =  0.8 


k 

C(pr) 

oc ,   . 
Cq  (pr) 

a°^(sec) 

0.8 

.9232 

.8680 

.8009 

0.9 

.8543 

.8108 

.6289 

1.0 

.8140 

.8140 

.5301 

Shape  Curve  for  9  =  45?9,  x°^  =  3.476, 
X^^  =  4.047,  X  =  0.8 


k 

-r<p-) 

C(p-> 

a°^{sec) 

0.6 

.9293 

.8541 

.6786 

.8313 

0.65 

.4568 

.9789 

0.7 

.7810 

.6877 

.3700 

0.8 

.6862 

.6092 

.2767 

0.9 

.5747 

.2245 

1.0 

.5828 

.1913 
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oc 
with  k(a   ,  q  )  taken  from  the  x  =  .8  tables  was  also 

derived  and  plotted  (Table  XIV  and  Figure  7) .   An  inspection 

of  this  plot  shows  that  for  0   =50°  and  0^  =  45?9,  the 
■^  e  e 

intersection  of  the  shape  curves  for  the  two  eclipses  in 
both  cases  lies  above  the  depth  curve.   This  implies  the 
possibility  of  a  solution  from  the  shapes  alone,  abandoning 
the  depth  curve  and  therefore  the  assumption  that  L,  +  L2=  1. 
This  so-called  third  light  solution  is  indeed  one  way  to 
produce  a  theoretical  light  curve  that  will  fit  the  observed 
curve.   There  is,  however,  no  real  justification  for  assuming 
the  presence  of  this  third  light  in  the  present  case  and 
therefore,  this  possibility  was  rejected. 

It  had  been  noticed  from  the  first  plotting  of  this 
rectified  light  curve  that  a  number  of  the  points  in  the 
primary  and  secondary  eclipses  were  similar  (i.e.  for  a  given 
value  of  0,  I"  was  nearly  the  same  for  both  eclipses) .   In 
attempting  solutions  for  k  =  1,  this  similarity  became 
even  more  apparent  and  it  was  decided  to  superimpose  the 
plots  of  the  two  eclipses  (see  Figure  8) .   The  result  was 
remarkable,  from  the  shoulders  down  to  a  depth  of  approxi- 
mately n  =  0.8,  the  two  eclipses  were  virtually  identical; 
deeper  than  this  point,  there  was  a  sharp  divergence.   Since 
this  type  of  behavior  is  not  possible  in  a  system  described 
by  the  Russell  model,  either  the  points  near  mid-eclipse  had 
to  be  completely  abandoned,  or  the  rectification  itself 
might  be  at  fault,  and  a  new  rectification  could  be  tried. 
It  was  decided  to  try  to  find  a  solution  on  the  other 
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TABLE  XIV 


oc 
Depth  Curve  for  x  =  0.8,  1-ii   =0.16, 


1   - 

■    il^^   =   0.18 
o 

^o 

o 

0.195 

1.00 

0.2 

.98 

0.25 

.82 

0.3 

.713 

0.4 

.58 

0.5 

.50 

0.6 

.447 

0.7 

.408 

0.8 

.38 

0.9 

.357 

1.0 

.34 

.425 
.433 
.510 
.565 
.650 
.724 
.788 
.848 
.901 
.952 
1.000 
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rectification  mentioned  above,  i.e.  the  one  whose  coefficients 
were  derived  by  the  least  squares  method  of  Fourier  analysis. 

Solutions  from  Least  Squares  Rectification 

There  are  several  differences  apparent  in  looking 
at  the  two  curves  produced  by  the  two  different  sets  of 
rectification  coefficients.   Outside  the  eclipses,  the 
residuals  are  more  or  less  evenly  distributed  throughout 
the  whole  curve  in  the  least  squares  rectification,  with 
the  sum  of  the  squares  of  the  residuals  smaller  than  that 
from  the  graphical  rectification,  as  expected.   Inside  the 
eclipses  the  two  curves  are  no  longer  virtually  identical 
from  the  shoulders  down  to  an  approximate  level  of  n  =  0.8; 
in  general,  for  any  given  phase  the  primary  is  deeper  than 
the  secondary.   The  greatest  change  in  the  new  rectification, 
however,  is  in  the  depths  of  the  two  eclipses.   Both 
eclipses  are  significantly  shallower,  the  depth  of  the 
primary  eclipse  going  from  about  0.18  to  0.16  and  the  depth 
of  the  secondary  eclipse  going  from  about  0.16  to  0.145. 

A  nomographic  solution  of  this  second  rectified 
light  curve  was  then  attempted.   The  light  curve  was  rec- 
tified in  phase  with  a  value  of  z  based  on  N  =  3.2  corre- 
sponding to  a  darkening  x  =  0.8  and  the  solution  was  tried 
on  the  X  =  0.8  nomograph.   The  value  of  x  =  0.8  was  chosen 
for  convenience  since  there  is  neither  a  nomograph  nor  a 
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set  of  X  tables  for  complete  darkening  and  the  previous 

explorations  of  the  light  curve  seemed  to  indicate  a  large 

value  for  the  darkening.   The  depths  of  the  eclipses 

were  chosen  as  1  -  l^^  =   0.16  and  1  -  ilj®^  =  0.14  and  the 

o  o 

following  values  were  read  off  the  plot  of  the  rectified 


light 

curve 

(Figure  4) : 

n 

eP^ 

sin^oP^ 

xP" 

gsec 

.  2p,sec 
sm  9 

^sec 

0.2 

30?0 

.25000 

2.038 

31?5 

.27300 

2.032 

0.5 

20.5 

.12265 

1 

21.5 

.13432 

1 

0.8 

12.0 

.043227 

.352 

13.5 

.054497 

.406 

Since  x    (n  =0.8)  >  x   (n  =  0.8),  it  seems  that  the 
primary  eclipse  is  a  transit  and  the  secondary  eclipse  is 
an  occultation  as  before.   The  values  needed  for  the  depth 
line  on  the  nomograph  are 


(1  -  1°"")    +    a   -    H^)    =  0.30 


1  -  t^^ 

=  0.1860 

£°^ 
o 


1  -  £^^    ,   1  -  ^"^ 
O   j^  1         o 


+  7^ =   0.1893 


jOC      50    pOC 
o  o 


Taking  these  values  and  the  x (n  =  0.8)  values  given  above 
to  the  X  =  0.8  nomograph,  an  intersection  of  the  depth  line 

and  both  x  contours  was  found  with  k  =  0.45,  p  =  -0.75. 

oc  tr 

a^  =  0.9395,  and  a^  =  0.8522.   Using  the  x  tables,  the 

following  points  on  the  light  curve  were  derived  from  these 


46 


parameters. 

n 

X?^ 

2 

sin  9 

G 

tr 

X 

0.0 
0.2 
0.5 
0.8 

3.037 
1.780 
1 
.407 

.44706 
.26202 
.14720 
.059912 

42?0 
30.8 
22.6 
14.2 

3.645 
1.967 
1 
.352 

sin^e     0 

.44706  42?0 

.24125  29.4 

.12265  20.5 

.043173  12.0 

The  situation  here  is  quite  similar  to  that  found  with 
the  previous  rectification.   The  computed  curve  is  defi- 
nitely narrower  at  the  shoulders  than  the  observed  curve. 
The  fit  of  primary  eclipse  from  the  half-way  point  down  is 
fairly  satisfactory.   The  lower  half  of  the  computed  curve 
for  the  secondary  eclipse  is  considerably  wider  than  the 
observed  curve. 

These  difficulties  are  basically  the  same  as  those 
encountered  in  the  solution  based  on  the  first  rectification. 
Furthermore,  the  situation  here  is  one  common  to  solutions 
of  many  systems  with  relatively  shallow  eclipses,  i.e. 
solutions  with  a  wide  range  of  values  of  k  differ  very 
little  in  the  light  curves  which  they  produce.   For  example, 
in  this  particular  case,  if  the  depths  of  the  two  eclipses 
and  the  value  of  the  external  tangency  point,  0  ,  are  fixed, 
then  solutions  for  the  cases:  1)  k  =  .65  primary-occultation, 
2)  k  =  .65  primary-transit,  3)  k  =  1  (see  Tables  XV,  XVI, 
and  XVII  and  Figure  4) ,  differ  in  the  whole  course  of  the 
light  curve  by  no  more  than  approximately  0?6.   This  means 
that  distinguishing  among  these  possible  solutions  is  quite 
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TABLE   XV 

Solution  for  x  =   0.8,    1   -    i^^  =   0.16,    1  -    5-^^^  =  0.145, 
pr-oc,   k  =  0.65,    p^=  -0.1625,    a°^=  0.55,a^^=^v4548 


n 

X°" 

sin^e 

0 

tr 

X 

2 

sin  0 

0 

0.0 

3.515 

.48255 

44?0 

3.696 

.48255 

44?0 

0.1 

2.418 

.331 

35.2 

2.507 

.32731 

34.9 

0.2 

1.935 

.26564 

31.0 

1.986 

.25929 

30.6 

0.4 

1.264 

.17353 

24.6 

1.275 

.16646 

24.1 

0.5 

1 

.13728 

21.7 

1 

.13056 

21.2 

0.6 

.764 

.104884 

18.9 

.758 

.098064 

18.3 

0.8 

.353 

.048461 

12.7 

.344 

.944913 

12.2 

0.9 

.171 

.023475 

8.8 

.165 

.021542 

8.4 

0.95 

.0840 

.011532 

6.2 

.0807 

.010536 

5.9 
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TABLE  XVI 


Solution  for  X  =  0.8,  1  -  l^^  =   0.16,  1  -  l^^^   =  0.145, 

'       o         '       o  ' 

pr-tr,  k  =  0.65,  p  =  -0.1271,  a°^  =  0.525,  a^^  =  0.4314 
^                  '^O  o  o 


n 

oc 

X 

sin  0 

0 

tr 

X 

sin  0 

0 

0.0 

3.505 

.48255 

44?0 

3.671 

.48255 

44?0 

0.1 

2.415 

.33249 

35.2 

2.497 

.32823 

35.0 

0.2 

1.934 

.26626 

31.1 

1.981 

.26040 

30.7 

0.4 

1.264 

.17402 

24.7 

1.274 

.16747 

24.1 

0.5 

1 

.13768 

21.8 

1 

.13145 

21.3 

0.6 

765 

.105321 

18.9 

.758 

.099638 

18.4 

0.8 

.354 

.048737 

12.8 

.345 

.045350 

12.3 

0.9 

.171 

.023542 

8.8 

.165 

.021689 

8.5 

0.95 

.0841 

.011578 

6.2 

.0811 

.010661 

5.9 
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TABLE  XVII 


Solution  for  X  =  0.8,  1  -  ii^^  =  0.16,  1  -  i^^^  =   0.145, 
k  =  1.0,  p^  =  0.1540,  a  =  0.305 


n  X 

0.0  3.626 

0.1  2.475 

0.2  1.969 

0.4  1.271 
0.5  1 

0.6  .760 

0.8  .347 

0.9  .167 

0.95  .0818 


2 

sin  0 

0 

.48255 

44?0 

.32938 

35.0 

.26204 

30.8 

.16915 

24.3 

.133081 

21.4 

.101142 

18.5 

.046179 

12.4 

.022225 

8.6 

.010886 

6.0 
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difficult. 

Because  of  the  above  difficulties,  a  provisional' 
solution  with  the  parameters  1  -  i^^  =   0.16,  1  -  £^®*^  =  0.145, 
k  =  1.0,  a^  =  .305,  and  p^  =  .1540  was  chosen.   The  computed 
curve  was  plotted  (Figure  4)  and  values  of  observed  minus 
computed  (0  -  c)  for  each  observed  point  were  derived  graph- 
ically.  These  (0  -  C)  values  are  listed  in  Table  XVIII  and 
plotted  as  a  function  of  phase  in  Figure  9 .   From  an  inspec- 
tion of  this  plot,  it  appears  that  a  higher  order  cosine 
term  (e.g.  cos  9  0)  is  present,  especially  in  the  primary 
eclipse. 

The  possibility  of  using  Kitamura's  method  for 
the  solution  of  eclipsing  binary  light  curves  (IT)  was  also 
explored.   Kitamura's  method  has  the  advantage  of  using  the 
whole  light  curve  to  find  a  provisional  solution  rather 
than  using  a  few  selected  points  as  is  done  in  the  nomo- 
graphic method.   Since  both  methods  of  solution  are  based 
upon  the  same  geometrical  model  for  the  eclipses,  their 
final  results  should  agree.   Even  though  the  two  methods 
might  produce  different  provisional  solutions,  a  careful 
and  thorough  analysis  based  upon  these  different  preliminary 
solutions  should  finally  produce  the  same  final  solution. 
Any  difference  in  the  results  has  to  be  caused  by  improper 
application  of  one  of  the  methods.   Since  solutions  using 
Kitamura's  method  should  not  differ  from  those  derived 
using  the  Russell-Merrill  method,  Kitamura's  method  was 
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TABLE  XVIII 
0-C's  from  Solution  for  k  =  1,  p.  =  0.1540,  a  =  0.305. 

0  O-CO)           O-C{180  +0) 

0.000  .007                .017 

2.945  -.003                .007 

5.885  -.007                .000 

8.816  -.003  -.002 

11.735  .003  -.002 

14.637  .000  -.004 

17.519  .003                .002 

20.377  .011                .007 

23.208  .006                .005 

26.011  .004                .002 

28.781  .001                .001 

34.221  -.007                .003 

39.518  -.006  -.001 
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not  pursued  any  further. 

Before  going  on  to  find  the  geometrical  elements 
corresponding  to  these  possible  solutions,  a  further 
exploration  of  the  possibility  that  the  eclipses  are  com- 
plete might  be  in  order.   By  making  this  assumption,  it 
is  possible  to  find  a  solution  from  the  shapes  of  the 
eclipses  alone  (i.e. ,  by  not  assuming  that  L,  +  L^  =  1) . 
While  this -introduction  of  third  light  seems  somewhat 
unwarranted,  it  is  presented  here  for  the  sake  of  com- 
pleteness.  The  shapes  of  the  eclipses  were  used  to  derive 
the  solution  given  in  Table  X IX  and  plotted  in  Figure  4 . 
This  solution  assumes  that  the  eclipses  are  complete  and, 
in  fact,  takes  the  limiting  situation  of  central  eclipses 
(i.e.,  i  =  90°).   These  assumptions  require  that  more  than 
half  of  the  total  light  of  the  system  come;  from  some 
unknown  third  body  (L_  =  0.588).   This  solution  seems  to 
fit  the  light  curve  about  as  well  as  the  previously  men- 
tioned solutions. 

In  rectifying  the  phase,  the  ellipticity  z  is 
a  function  of  a  parameter  N  which  goes  as 

N  =  (15  +  X)  (1  +  y) 
15  -  5x 

where  x  is  the  limb  darkening  and  y  is  the  gravity  darkening 
of  the  star  being  eclipsed.   In  the  rectification  of  the 
phase  described  above,  a  value  of  N  =  2.6  corresponding  to 
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TABLE   Xix: 

Solution  for  X  =  0.8,    1   -    Ji^^  =   0.165,    1   -    ii^®^  =  0.145, 
k   =   0.70,    p^   =   -1.429,    a°^   =   1.0,    a^^  =   1.067, 


T  =  0.558590,  L^ 

=  0.588. 

n 

oc 

X 

2 

sin  0 

0 

tr 

X 

.  2 

sm  0 

0 

0.0 

3.339 

.48255 

44?0 

4.481 

.48255 

44?0 

0.1 

2.281 

.32965 

35.0 

2.833 

.30508 

33.5 

0.2 

1.834 

.26505 

31.0 

2.163 

.23293 

28.9 

0.4 

1.231 

.17790 

24.9 

1.307 

.14075 

22.0 

0.5 

1 

.14452 

22.3 

1 

.107688 

19.2 

0.6 

.798 

.11533 

19.9 

.745 

.080228 

16.5 

0.8 

.452 

.065323 

14.8 

.350 

.037691 

11.2 

0.9 

.297 

.042922 

12.0 

.196 

.021107 

8.4 

0.95 

.217 

.031361 

10.2 

.1142 

.012298 

6.4 

1.0 

.1040 

.015030 

7.0 

0.937=n, 

.1396 

.015033 

7.0 
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a  limb  darkening  of  x  =  0.6  implies  that  a  gravity  darkening 
of  y  =  1.0  was  chosen.   Dr.  S.  Rucinski  (6)  has  suggested 
that  a  value  of  gravity  darkening  y  =  0.32  might  be  more 
appropriate  for  this  system.   This  possibility  was  explored 
and  it  was  found  to  have  no  significant  effect  on  the 
problems  encountered  in  finding  a  solution  for  this  system. 
Therefore,  it  was  decided  to  proceed  with  the  solutions 
already  derived. 

Orbital  Elements 

The  next  step  was  to  derive  orbital  elements  for 
the  eclipsing  system  from  the  geometrical  eclipse  parameters. 
These  elements  can  most  easily  be  derived  from  the  formulas 
given  by  Merrill  (5)  : 


2 

o  sm  0 

a^  =  ^ 


(1  +  k)^  -  (1  +  kp  )^cos^0 
o       e 


a  =  ka 
s     g 


cos  i^  =  (1  +  kp^)ag 


T      ,     „OC  /  oc 
L  =1-1      /a 
s        o  '  o 


L^  =  1  -  L 
g       s 
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The  derived  orbital  elements  for  the  various  solutions  are 
listed  in  Table  XX. 


Mean  errors  (y  =  \/z(0  -  C)^/(n  -  1))  were  derived 
for  the  "observed"  points  from  each  of  the  four  possible 
solutions.   These  mean  errors  are  listed  in  Table  XX.   The 
mean  error  has  almost  the  same  value  for  all  three  of  the 
partial  eclipse  solutions  (with  no  third  light  added) .   The 
value  of  the  mean  error  for  the  complete  eclipse  solution 

(with  third  light)  is  slightly  smaller  than  the  values  of 
the  mean  error  for  the  other  three  solutions.   Since  these 
four  solutions  with  very  different  values  of  the  geometrical 
elements  produce  computed  light  curves  which  fit  the  ob- 
served light  curve  more  or  less  equally  well,  it  would  not 
be  meaningful  to  apply  the  method  of  differential  corrections 
to  find  a  least  squares  solution  of  the  light  curve.   In 
the  absence  of  a  formal  least  squares  solution,  it  is  not 
possible  to  determine  mean  errors  for  each  of  the  geome- 
trical elements.   From  the  results  presented  in  Table  XX, 
however,  it  is  possible  to  estimate  a  reasonable  range  of 
each  element.   For  the  three  partial  eclipse  solutions 

(without  third  light)  the  approximate  ranges  of  these 

elements  are:  .38  <ja   <  .46;  .29  <  a^  <  .38;  64°  <  i^  <  66°; 

g  s  r 

153  <  L  <  .72;   and  .27  <  L   <  .48.    With  the  addition 
g  s 

of  third  light  the  inclination  can  increase  up  to  the 
limiting  value  of  90°  and  solutions  are  possible  anywhere 
within  this  range  of  inclinations.   In  the  solutions 
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TABLE  XX 
Orbital  Elements  for  Solutions  of  First  Synthetic  Light  Curve 


pr-oc 

pr-tr 

k 

.65 

110: 

.65 

.70 

o 

.5500 

.3050 

.5250 

1 

.0 

tr 
^'o 

.4548 

.3050 

.4314 

1 

.067 

Po 

-.1625 

+.1540 

-.1271 

-1 

.429 

a 

g 

.457 

.382 

.459 

.409 

^s 

.297 

.382 

.299 

.286 

cos  i^ 

.40890 

.44083 

.42139 

0 

i 
r 

65?9 

63?9 

65?1 

90° 

L 

g 

.709 

.525 

.724 

.277 

^s 

.291 

.475 

.276 

.135 

h 

0 

0 

0 

.588 

y 

.0045 

.0044 

.0044 

.0037 
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presented  here,  the  value  of  0  was  chosen  as  44°  in  order 

e 

to  best  fit  the  whole  light  curve.   The  observations  on 

the  shoulders  of  the  eclipses  alone  indicate  that  this 

value  might  be  increased  by  several  degrees.   The  effect 

of  increasing  0   is  to  increase  the  values  of  a  and 
e  g 

a  (i.e.,  to  increase  the  sizes  of  the  stars  relative  to 
their  separation)  and  to  decrease  the  inclination  of  the 
orbit. 


Figures  of  the  Components 

The  figures  of  the  components  can  now  be  computed 
from  the  formulas  given  by  Merrill  (10) : 


r  +  (0.17  +  1.19  m  /m  )r^   =   a 
g  s^  g'  g    g 


bg=  ag  -  1.53(m^/mg)r^ 


c^  =  3r  -  a  -  b 

g    g   g   g 


These  formulas  require  that  the  mass  ratio  of  the  compo-  ■ 
nents  be  known.   Since  this  mass  ratio  was  not  known,  a 
more  indirect  method  was  used  to  set  limits  on  both  the 
figures  of  the  components  and  their  mass  ratios.   This 
method  assumes  that  no  dimension  of  the  star  may  exceed  the 
size  of  the  Roche  lobe  it  is  contained  within.   The  dimen- 
sions of  the  Roche  lobes  for  various  mass  ratios  are 
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tabulated  in  the  literature  (e.g.,  Kopal  (12)).   The 
figures  of  the  components  were  computed  for  a  number  of 
different  mass  ratios  for  each  of  the  solutions.   (The 
two  solutions  for  k  =  .65  were  combined  for  this  purpose 
since  the  values  of  a  and  a  in  each  were  nearly  equals ) 
These  figures  were  compared  with  the  size  of  the  Roche  lobe 
in  order  to  determine  the  limiting  mass  ratios  which  would 
satisfy  the  criterion  that  the  sizes  of  the  stars  must  not 
exceed  the  sizes  of  their  Roche  lobes.   These  limiting 
mass  ratios  and  the  figures  of  the  components  derived 
using  them  for  each  of  the  solutions  are  listed  in  Table 
XXI. 

Comparison  of  Solutions  with  Input  Parameters 

It  is:  now  possible  to  compare  the  solutions 

derived  above  with  the  parameters  used  to  generate  the 

light  curve.   The  light  curve  was  produced  using  Lucy's 

model  for  the  light  curves  of  W  Ursae  Majoris  stars.   In 

this  model  the  surfaces  of  the  two  stars  are  Roche  equipo- 

tential  surfaces  and  for  the  Rucinski  light  curve  the 

stars  share  a  common  envelope  which  fills  the  outer 

contact  surface.   The  constant  which  defined  this  surface 

is  C  =  3.5591  with  the  mass  ratio  q  =  m  /m  =  0.4.   The 

5   y 

Other  parameters  used  in  deriving  the  light  curve  are  the 
inclination  of  the  orbit  i  =  82°,  the  effective  temperature 
of  the  surface  T  =  5700°,  the  wavelength  of  the  light 
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TABLE  XXI 


Figures  of  the  Components 


k 

.65 

.65 

1.0 

1.0 

.70 

.70 

m^/mg 

.28 

.57 

.87 

1.15 

.25 

.96 

a 

g 

.458 

.458 

.382 

.382 

.409 

.409 

^g 

.439 

.429 

.361 

.357 

.397 

.381 

b 

g 

.442 

.428 

.359 

.354 

.400 

.378 

c 

g 

.417 

.401 

.342 

.335 

.382 

.356 

^s 

.298 

.298 

.382 

.382 

.286 

.286 

^s 

.273 

.283 

.357 

.361 

.263 

.278 

^s 

.267 

.281 

.354 

.359 

.257 

.276 

c„ 

.254 

.270 

.335 

.342 

.246 

.272 
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X  =  5500  A,  the  limb  darkening  u  =  0.6,  and  the  gravity 

darkening  y  =  0.32.   The  size  of  the  equipotential  surface 

(and  therefore  the  stars)  was  also  given.   These  values 

are  b  =  0.531,  c  =  0.482,  b  =0.363,  and  c  =0.337  :(60 
g  g  s  s 

These  parameters  can  be  compared  with  the  results 
of  the  solutions  as  listed  in  Table  XX  and  Table  XXI.   The 
inclination  of  the  orbit  suggests  that  a  solution  using 
the  Russell-Merrill  method  which  has  this  inclination  will 
include  the  presence  of  some  third  light.   The  figures  of 
the  components  (with  q  =  0.4)  for  such  a  solution  can  be 
of  roughly  the  same  shape  in  the  yz-plane  (the  plane  con- 
taining the  b  and  c  axes  of  the  components)  as  the  shapes 
of  the  lobes  of  the  outer  contact  surface  in  the  yz-plane. 
The  Russell-Merrill  dimensions  of  the  components  in  this 
plane  relative  to  the  separation  of  the  components  are 
considerably  smaller  than  the  dimensions  of  the  lobe  of 
the  outer  contact  surface  in  the  yz-plane.   Although  it  is 
possible,  as  previously  mentioned,  to  increase  slightly 
the  dimensions  of  the  components  by  increasing  the  values 
of  0   in  the  solution,  it  is  not  possible  to  increase  the 
dimensions  of  the  components  to  anywhere  near  the  size  of 
the  outer  contact  surface  and  still  have  the  value  of  0 
compatible  with  the  light  curve.   Using  the  Russell  model 
of  similar  prolate  ellipsoids,  it  is  possible  to  propose 
a  system  consisting  of  two  components  whose  dimensions  in 
the  yz-plane  are  nearly  as  large  as  the  dimensions  of  the 
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lobes  of  the  outer  contact  surface  in  the  yz-plane.   These 
prolate  ellipsoids  would  be  in  contact  (i.e.,  a  +  a^  =  1) . 
There  are  a  number  of  objections  to  this  proposed  system. 
First  of  all,  it  is  not  possible  for  real  stars  in  contact 
to  have  ellipsoidal  shapes.   Even  if  such  a  system  were 
physically  possible,  it  would  produce  eclipses  with  e^  =  90°, 
in  contrast  with  the  values  (44°<  Q^  <   50°)  given  by  the 
light  curve.   Even  the  physically  realistic  system  with 
the  two  components  filling  their  Roche  lobes  (inner  contact 
surfaces)  produces  eclipses  that  are  wider  (for  i  =  82°) 
than  those  in  the  light  curve  (12).   In  this  case  the  sizes 
of  the  components  in  the  yz-plane  are  also  considerably  smaller 
than  the  sizes  of  the  lobes  of  the  outer  contact  surface  in 
the  yz-plane. 

As  the  above  discussion  indicates,  the  solutions 
produced  by  applying  the  Russell-Merrill  method  to  the  light 
curve  give  very  different  values  for  the  orbital  elements 
than  the  parameters  used  to  generate  the  light  curve.   This 
is  especially  true  if  third  light  is  not  included  in  the 
solution.   Then,  not  only  are  the  stars  much  smaller  than 
the  outer  contact  surface,  but  the  inclination  of  the  orbit 
is  also  quite  different  from  the  inclination  used  in  de- 
riving the  light  curve.   Leaving  aside  the  question  of 
third  light,  since  the  amount  of  third  light  cannot  be 
determined  from  the  light  curve  alone,  it  is  apparent  that 
the  Russell-Merrill  solution  of  the  light  curve  is  com- 
pletely unsuccessful  in  predicting  the  elements  used  to 
derive  the  light  curve  from  Lucy's  model.   Thus  these  two 
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models  for  the  light  curves  of  W  Ursae  Ma j oris  stars  seem 
to  be  incompatible.   While  the  models  can  produce  similar 
light  curves,  both  of  which  resemble  the  light  curves 
observed  for  W  Ursae  Majoris  type  stars,  the  orbital  elements 
which  they  use  to  produce  these  similar  light  curves  are 
quite  different.   Since  most  of  the  solutions  of  the  light 
curves  of  W  Ursae  Majoris  type  stars  which  are  in  the 
literature  are  based  upon  the  geometrical  model  lying 
behind  the  Russell-Merrill  method  of  solution  (or,  equiv- 
alently,  other  methods  based  upon  the  same  model) ,  much 
caution  is  in  order  in  looking  at  these  solutions.   If 
Lucy's  model  represents  the  actual  physical  situation  for 
a  W  Ursae  Majoris  type  star,  then  most  of  the  solutions  in 
the  literature  probably  do  not  realistically  represent  the 
actual  physical  situation  present  in  the  W  Ursae  Majoris 
type  stars. 


CHAPTER  III 
SECOND  SYNTHETIC  LIGHT  CURVE 

A  second  synthetic  light  curve  for  a  W  Ursae  Ma j oris 
type  eclipsing  binary  was  generated  from  a  theoretical  model 
by  Dr.  R.  E.  Wilson  (9) .   The  "observational"  data  for  this 
light  curve  consisted  of  51  values  of  the  differential 
stellar  magnitude  as  a  function  of  phase,  with  the  phases 
given  in  intervals  of  one-hundredth  of  the  period  from 
phase  0.00  to  phase  0.50.   The  remaining  half  of  the  light 
curve  from  phase  0.50  to  phase  1.00  was  assumed  to  be 
symmetrical.   This  "observational"  data  is  plotted  in  Figure 
10  and  listed  in  Table  XXII.   (A  normalization  constant 
equal  to  0.2049  was  subtracted  from  each  of  the  magnitudes;) 
Also  listed  in  this  table  are  the  corresponding  light  values 
and  phase  angles  in  degrees  for  each  data  point.   The  syn- 
thetic light  curve  was  to  be  treated  as  observational  data 
and  a  solution  was  to  be  obtained  by  using  the  standard 
Russell-Merrill  method  of  solution  for  light  curves  of 
eclipsing  binaries.   No  additional  information  about  the 
nature  of  the  derivation  of  this  light  curve  was  to  be  used 
in  the  solution. 
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TABLE  XXII 
Observations  of  Second  Synthetic  Light  Curve 


PHASE 

Ain-0.2049 

I 

e 

0.00 

0.6900 

0.5297 

0.0 

0.01 

0.6893 

0.5300 

3.6 

0.02 

0.6857 

0.5318 

7.2 

0.03 

0.6398 

0.5547 

10.8 

0.04 

0.5725 

0.5902 

14.4 

0.05 

0.5012 

0.6303 

18.0 

0.06 

0.4333 

0.6709 

21.6 

0.07 

0.3718 

0.7100 

25.2 

0.08 

0.3176 

0.7464 

28.8 

0.09 

0.2707 

0.7793 

32.4 

0.10 

0.2308 

0.8085 

36.0 

0.11 

0.1972 

0.8339 

39.6 

0.12 

0.1691 

0.8558 

43.2 

0.13 

0.1457 

0.8744 

46.8 

0.14 

0.1252 

0.8911 

50.4 

0.15 

0.1060 

0.9070 

54.0 

0.16 

0.0881 

0.9221 

57.6 

0.17 

0.0712 

0.9365 

61.2 

0.18 

0.0557 

0.9500 

64.8 

0.19 

0.0419 

0.9621 

68.4 

0.20 

0.0298 

0.9729 

72.0 

0.21 

0.0197 

0.9820 

75.6 

0.22 

0.0115 

0.9895 

79.2 

0.23 

0.0057 

0.9948 

82.8 

0.24 

0.0016 

0.9985 

86.4 

0.25 

0.0000 

1.0000 

90.0 

0.26 

0.0006 

0.9994 

93.6 

0.27 

0.0037 

0.9966 

97.2 

0.28 

0.0087 

0.9920 

100.8 

0.29 

0.0160 

0.9854 

104.4 

0.30 

0.0256 

0.9767 

108.0 

0.31 

0.0374 

0.9661 

111.6 

0.32 

0.0510 

0.9541 

115.2 

Table  XXII      continued. 
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PHASE 

Am-0.2049 

I 

e 

0.33 

0.0667 

0.9404 

118.8 

0.34 

0.0841 

0.9255 

122.4 

0.35 

0.1028 

0.9097 

126.0 

0.36 

0.1231 

0.8928 

129.6 

0.37 

0.1459 

0.8743 

133.2 

0.38 

0.1715 

0.8539 

136.8 

0.39 

0.1998 

0.8319 

140.4 

0.40 

0.2322 

0.8075 

144.0 

0.41 

0.0694 

0.7803 

147.6 

0.42 

0.3115 

0.7506 

151.2 

0.43 

0.3621 

0.7164 

154.8 

0.44 

0.4172 

0.6810 

158.4 

0.45 

0.4779 

0.6439 

162.0 

0.46 

0.5426 

0.6067 

165.6 

0.47 

0.6066 

0.5720 

169.2 

0.48 

0.6574 

0.5458 

172.8 

0.49 

0.6672 

0.5409 

176.4 

0.50 

0.6693 

0.5399 

180.0 
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Rectification 

As  in  the  previous  synthetic  light  curve  solution, 
the  first  step  is  an  analysis  of  the  light  outside  the 
eclipses  in  order  to  rectify  the  light  curve  to  the  spher- 
ical model.   This  analysis  was  done  by  a  least-squares 
Fourier  analysis  of  the  light  outside  of  eclipse  (48°  to 
132°) .   The  results  of  this  analysis  are  listed  in  Table 
XXIII.   The  light  values  were  then  rectified  by  the 
formula : 


I+C  +Ct  cos  e+C^  cos  2e-A-  cos  3e-A,  cos  49 
jii  _    o   1 2 3 4 

A  +C  +  (A^+C^)  cos  26 
O   O     2   2 


and  the  phase  angles  were  rectified  by  the  formula 


2         2  2 

sin  0  =  sin  6/(1  -  z  cos  6) 


where  z  =  -4 (A2  -  C2)/N(Aq  -  C^  -  A2  +  C2)  with  N  =  2.2, 
corresponding  to  a  limb  darkening  x  =  0.4.   The  rectification 
coefficients  for  reflection  were  determined  from  the  statis- 
tical formula  given  by  Russell  and  Merrill  (2) :  C  =  0.072 

sin^e      =   0.0398,    C,    =   -A,,    and   C-   =   0.024    sin   6      =   0.0133, 
e  11  2  e 

with  6   =48°.   Following  Merrill  (10),  the  A^  and  A^ 

terms  were  removed  by  subtraction.   The  rectified  points  are 

listed  in  Table  XXIV  and  plotted  in  Figure  11. 
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TABLE  ;  XXI 1 1 

Rectification  Coefficients  for  Second 
Synthetic  Light  Curve 


A' 
o 

=   0.86563 

+ 

0.00012 

^1 

=   0.00390 

+ 

0.00009 

^2 

=  -0.13528 

+ 

0.00018 

^3 

=   0.00382 

+ 

0.00005 

^4 

=  -0.00098 

+ 

0.00007 

OBSERVED 

THEORETICAL 

o-c 

.,  e 

INTENSITY 

INTENSITY 

0.8911 

0.8910 

0.0001 

50.4000 

0.9070 

0.9069 

0.0001 

?;54.'0000 

0.9221 

0.9222 

-0.0001 

57.6000 

0.9365 

0.9366 

-0.0001 

61.2000 

0.9500 

0.9500 

-0.0000 

64.8000 

0.9621 

0.9622 

-0.0001 

68.4000 

0.9729 

0.9729 

0.0000 

72.0000 

0.9820 

0.9820 

-0.0000 

75.6000 

0.9895 

0.9894 

0.0001 

79.2000 

0.9948 

0.9949 

-0.0001 

82.8000 

0.9985 

0.9984 

0.0001 

86.4000 

1.0000 

0.9999 

0.0001 

90.0000 

0.9994 

0.9994 

0.0000 

93.6000 

0.9966 

0.9967 

-0.0001 

97.2000 

0.9920 

0.9920 

-0.0000 

100.8000 

0.9854 

0.9853 

0.0001 

104.4000 

0.9767 

0.9767 

0.0000 

108.0000 

0.9661 

0.9662 

-0.0001 

111.6000 

0.9541 

0.9541 

0.0000 

115.2000 

0.9404 

0.9405 

-0.0001 

118.8000 

0.9255 

0.9256 

-0.0001 

122.4000 

0.9097 

0.9096 

0.0001 

126.0000 

0.8928 

0.8928 

0.0000 

129.6000 
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TABLE  XXIV 
Rectified  Second  Synthetic  Light  Curve 

I  I"  0 


0.5297 

0.7353 

0.000 

0.5300 

0.7347 

4.233 

0.5318 

0.7341 

8.453 

0.5547 

0.7583 

12.647 

0.5902 

0.7960 

16.805 

0.6303 

0.8370 

20.917 

0.6709 

0.8760 

24.972 

0.7100 

0.9106 

28.965 

0.7464 

0.9396 

32.889 

0.7793 

0.9624 

36.741 

0.8085 

0.9792 

40.518 

0.8339 

0.9905 

44.219 

0.8558 

0.9970 

47.845 

0.8744 

0.9996 

51.398 

0.8911 

1.0001 

54.881 

0.9070 

1.0002 

58.298 

0.9221 

1.0000 

61.653 

0.9365 

0.9999 

64.950 

0.9500 

1.0000 

68.196 

0.9621 

1.0000 

71.397 

0.9729 

1.0001 

74.559 

0.9820 

1.0000 

77.687 

0.9895 

1.0002 

80.789 

0.9948 

1.0000 

83.870 

0.9985 

1.0001 

86.939 

1.0000 

1.0001 

90.000 

0.9994 

1.0001 

93.062 

0.9966 

0.9999 

96.130 

0.9920 

1.0000 

99.212 

0.9854 

1.0002 

102.314 

0.9767 

1.0001 

105.442 

0.9661 

0.9999 

108.604 

0.9541 

1.0001 

111.804 
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Table  XXIV  continued, 


I" 


0.9404 

1.0000 

115.051 

0.9255 

1.0000 

118.348 

0.9097 

1.0002 

121.703 

0.8928 

1.0001 

125.119 

0.9743 

0.9989 

128.602 

0.8539 

0.9958 

132.156 

0.8319 

0.9908 

135.782 

0.8075 

0.9826 

139.483 

0.7803 

0.9701 

143.260 

0.7506 

0.9533 

147.112 

0.7164 

0.9291 

151.036 

0.6810 

0.9012 

155.028 

0.6439 

0.8686 

159.084 

0.6067 

0.8331 

163.195 

0.5720 

0.7979 

167.354 

0.5458 

0.7708 

171.548 

0.5409 

0.7681 

175.768 

0.5399 

0.7681 

180.000 
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Solutions 

An  examination  of  the  rectified  light  curve  reveals 
several  things.   The  most  obvious  is  that  the  eclipses  are 
apparently  complete.   Furthermore  it  appears  that  the 
primary  eclipse  is  the  occultation  (total)  eclipse  and  the 
secondary  eclipse  is  the  transit  (annular)  eclipse.   Closer 
examination  reveals  a  very  slight  brightening  in  the  center 
of  the  primary  eclipse.   The  flatness  of  both  eclipses 
indicates  that  the  limb  darkening  is  probably  small.   The 

depths  of  the  two  eclipses  were  chosen  as  1  -  5,^  =  0.265 

Sep 
and  1  -  a  =   0.233  for  the  purpose  of  obtaining  a  pre- 

liminary nomographic  solution.   Then  the  following  values 
were  read  off  the  plot  of  the  rectified  light  curve: 

0P^      sin^e     xP""     e^^^     sin^Q   f^"" 

0.2     33.5    •=^..30463     1.841     33.1     .29823    1.986 
0.5     24.0     .16543      1      22.8     .15017      1 
0.8     15.8     .074137     .448     14.4     .061847    .412 

X^^(n  =  0.8)  >  x^^*^  (n  =  0.8),  confiinning  that  the  primary 
eclipse  is  the  occultation  eclipse  and  the  secondary 
eclipse  is  the  transit  eclipse.   These  x  values  along  with 

a  depth  line  1  -  i^^/SL^^^  =   0.3176  were  taken  to  the  total 

^  o   o 

eclipse  portion  of  the  x  =  0.4  nomographs  and  a  satisfactory 

intersection  of  the  three  curves  was  found.   This  intersection 

oc 
implied  a  solution  with  the  values  k  =  0.545,  a   =1.0, 

a^^=  1.016,  p  =  -1.10,  and  t=  0.3126.   The  results  of  this 
o         ^o 

solution  are  listed  in  Table  XXV  and  plotted  in  Figure  11. 
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TABLE  XXV 


Solution  for  x  =  0.4,  k  =  0.545,  a °^  =  1.0,  a ^^  =  1.016, 


=  -1.10,  T=  0.3126,  1  -  ll"^  =   0.265,  1  "  i?-|  =  0 


=  0.233 


n 

X°" 

sin  0 

9 

x" 

•  2^ 
sm  0 

0 

0.0 

2 

.977 

.50000 

45?0 

3.364 

.50000 

45?0 

0.1 

2 

.180 

.36613 

37.2 

2.395 

.35597 

36.6 

0.2 

1 

.785 

.29979 

33.2 

1.921 

.28552 

32.3 

0.4 

1 

.224 

.20557 

27.0 

1.256 

.18668 

25.6 

0.5 

1 

.16795 

24.2 

1 

.14863 

22.7 

0.6 

.800 

.13436 

21.5 

777 

.11549 

19.9 

0.8 

.447 

.075074 

15.9 

.398 

.059155 

14.1 

0.9 

.279 

.046858 

12.5 

.231 

.034334 

10.7 

0.95 

.190 

.031911 

10.3 

.1426 

.021195 

8.4 

1.0 

.0623 

.010463 

5.9 

0.984= 

=n. 

.069 

.010255 

5.8 

la 


A  brief  inspection  of  this  plot  shows  that  the  fit 
is  reasohably  satisfactory  except  at  the  shoulders  of  the 
eclipses.   At  the  shoulders,  the  computed  curve  is  sig- 
nificantly narrower  than  the  "observed"  points.   This ,  in 
essence,  is  the  same  situation  that  was  found  in  the  solu- 
tion of  the  first  synthetic  light  curve.   Further  refine- 
ments of  the  solution,  consistent  with  the  assumption  of 
complete  eclipses  and  the  retention  of  the  depth  line 

-2—  =  0.3176,  will  not  significantly  improve  this 

*o 
situation.   It  seems  that  there  is  probably  some  incom- 
patibility between  the  Russell  model  of  prolate  spheroids  as 
rectified  to  the  spherical  model  and  the  theoretical  model 
which  was  used  to  generate  the  light  curves.   If  the  assump- 
tion that  the  theoretical  model  more  closely  approximates  p 
physical  reality  for  W  Ursae  Majoris  stars  is  warranted,  then 
the  Russell  model  has  serious  faults  when  applied  to  these 
systems.   It  is  well  known,  of  course,  that  the  Russell  model 
is  only  a  rough  approximation  to  the  very  close  binary  systems. 
The  real  question,  which  this  study  will  attempt  to  answer 
partially,  is  whether  a  solution  of  such  a  system  derived 
from  the  Russell  model  will  provide  an  adequate  representation 
of  the  true  orbital  elements  of  the  binary  system. 

One  way  to  somewhat  improve  the  fit  of  a  solution 
to  the' "observed"  points  is  to  abandon  the  depth  line  and 
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do  a  solution  only  from  the  shapes  of  the  eclipses.   This 
may  mean  the  introduction  of  the  light  of  a  third  body  to 
the  system.   In  the  present  case,  the  introduction  of  third 
light  (L-.  =  0.234)  yields  the  central  eclipse  solution  listed 
in  Table  XXVI  and  plotted  in  Figure  11.   This  solution  im- 
proves the  fit  somewhat,  though  the  shoulders  of  the  com- 
puted curve  are  still  too  narrow.   There  seems  to  be  no  sound 
observational  reason  for  adding  this  third  light.   In  about 
20%  of  the  solutions  of  eclipsing  binary  light  curves  in  the 
literature  (13) ,  some  third  light  is  present.   This  third 
light  solution  is  presented  here  in  order  to  demonstrate 
how  the  addition  of  third  light  affects  the  solution.   Unless 
there  is  some  real  physical  source  of  this  third  light,  how- 
ever, the  necessity  for  its  introduction  to  improve  the  so- 
lution only  serves  to  point  out  the  inadequacies  of  the 
model  on  which  the  solution  is  based. 

Orbital  Elements  and  Figures  of  the  Components 

These  two  solutions,  from  the  formulas  given  above, 
were  then  used  to  derive  orbital  parameters  and  limiting 
mass  ratios  for  the  assumption  that  neither  component  exceeds 
the  size  of  its  Roche  lobe.   The  results  of  these  derivations 
are  listed  in  Tables  XXVII  and  XXVIII. 

Comparison  of  Solutions  with  Input  Parameters 

It  is  now  possible  to  compare  the  solutions  derived 
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TABLE  XXVI 
Solution  for  X  =  0.4,  k  =  0.65,  a°^  =  1.0,  a^^  =  1.039, 


P 


=  -1.538,T  =  0.447537,  1   -    i^^  =   0.265> 


sec 


1  -  i^^^  =   0.233,  L3  =  0.234 


n 

oc 

X 

sin^e 

0 

tr 

X 

2 

sin  Q 

0 

0.0 

3.131 

.52965 

46?7 

3.600 

.52965 

46?7 

0.1 

2.250 

.38061 

38.1 

2.502 

.36812 

37.4 

0.2 

1.822 

.30821 

33.7 

1.976 

.29073 

32.4 

0.4 

1.229 

.20790 

27.1 

1.265 

.18612 

25.6 

0.5 

1 

.16916 

24.3 

1 

.14713 

22.6 

0.6 

.800 

.13533 

21.6 

.775 

.11403 

19.8 

0.8 

.462 

.078152 

16.2 

.413 

.060765 

14.3 

0.9 

.314 

.053116 

13.3 

.263 

.038695 

11.3 

0.95 

.240 

.040598 

11.6 

.1867 

.027469 

9.5 

1.0 

.141 

.023852 

8.9 

0.962= 

=^i 

.162 

.023835 

8.9 
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TABLE  XXVII 

Orbital  Elements  for  Solutions  of  Second 
Synthetic  Light  Curve 


k 

0.545 

0.65 

oc 
°'o 

1.0 

1.0 

tr 
^'o 

1.016 

1.039 

Po 

-1.10 

-1.538 

\ 

0.466:fi 

0.441 

a 
s 

0.254 

0.287 

cos  i 
r 

0.18663 

0 

^r 

79J2 

90° 

^g 

0.735 

0.501 

L 
s 

0.265 

0.265 

L. 

0 

0.234 

TABLE  XXVIII 
Figures  of  the  Components 


m^/nig 

.52 

.17 

.66 

.25 

a 

g 

.466 

.466 

.441 

.441 

r 

g 

.437 

.451 

.413 

.426 

b 

g 

.437 

.455 

.412 

.429 

c 

g 

.408 

.432 

.386 

.408 

a 
s 

.254 

.254 

.287 

.287 

r 
s 

.245 

.233 

.276 

.264 

■^s 

.244 

.228 

.274 

.258 

c 
s 

.237 

.217 

.267 

.247 
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above  with  the  parameters  used  to  generate  the  light  curve. 
The  theoretical  light  curve  is  a  trial  and  error  match  to 
Broglia's  observations  of  RZ  Comae  (14).   This  light  curve 
was  generated  by  Dr.  R.  E.  Wilson  from  the  parameters 
listed  in  Table  XXIX.    In  this  model,  as  in  the  model 
used  to  generate  the  first  synthetic  light  curve,  the  sizes 
of  both  the  components  exceed  the  sizes  of  their  Roche  lobes. 
In  this  case  the  boundary  of  the  components  lies  along  a 
common  envelope  which  is.  smaller  than  the  outer  contact 
surface.   Broglia's  observations  of  RZ  Comae  were  solved 
by  Binnendijk  (15) .   The  solution  of  the  synthetic  light 
curve  (without  third  light)  as  listed  in  Table  XXVII  and 
Table  XXVIII  is  fairly  close  to  Binnendijk 's  solution  of  the 
real  :  -  observations  of  RZ  Comae.   This  would  seem  to 
indicate  that  this  solution  is  the  one  which  follows  from 
the  Russell  model  and  that  it  is  not  greatly  sensitive  to 
effects  of  judgment  and  details  of  procedure.   However 
the  true  parameters  used  to  derive  the  synthetic  light 
curve  are  not  at  all  close  to  the  solution  given  above  or 
to  Binnendijk' s.   In  fact,  it  is  possible  to  derive  a 
theoretical  light  curve  from  the  true  parameters  using  the 
Russell  model.   Such  a  light  curve  would  look  very  different 
from  the  synthetic  light  curve  derived  from  the  same  para- 
meters using  Wilson's  model.   The  most  obvious  difference 
is   that  the  Russell  model  theoretical  light  curve  would 
have  significantly  wider  and  deeper  eclipses  than  those 
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TABLE  XXIX 


Input  Parameters  for  Second  Synthetic  Light  Curve 


i 

86?00 

^1 

0.3149 

J- 

^2 

0.6851 

^1 

0.40 

^2 

0.40 

r^^  (pole) 

0.2992 

r ■ (side) 

.3132 

r, (back) 
Xeff 

.3505 
5500  A 

^1 

1.00 

^2 

1.00 

^1 

5500K  (polar) 

^^2 

5563K(polar) 

"^2^^! 

2.200 

r^  (pole) 

0.4287 

r2  (side) 

.  .'4579 

k 

0.69 

r^  (back) 

.4880 

% 

5.449 

^2 

5.449 

^1 

1.00 

H 

1.00 

H 

0.000 
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present  in  the  given  synthetic  light  curve.   The  difference 
in  depth  is  especially  obvious  since  the  light  lost  in  the 
total  eclipse  is  equal  to  the  light  of  the  smaller  star.   In 
addition  the  annular  eclipse  in  the  theoretical  Russell 
Model  light  curve  derived  from  the  true  parameters  would  be 
deeper  than  the  total  eclipse,  the  reverse  of  the  situation 
in  the  synthetic  light  curve. 

A  comparison  of  the  solutions  derived  above  using 
the  Russell-Merrill  method  (Table  XXVII  and  Table  XXVIII) 
with  the  true  parameters  used  to  derive  the  synthetic  light 
curve  (Table  XX IX))  reveals  that  the  Russell-Merrill 
solution  of  the  light  curve  is  completely  unsuccessful  in 
predicting  the  elements  used  to  derive  the  light  curve 
from  Wilson's  model.   The  situation  here  is  similar  to  that 
found  in  the  solution  of  the  first  synthetic  light  curve. 
Wilson's  model,  like  Lucy's  model,  seems  to  be  incompatible 
with  the  Russell  model  as  applied  to  the  light  curves  of 
W  Ursae  Majoris  stars.   Again,  although  Wilson's  model  and 
the  Russell  model  will  produce  similar  light  curves,  both 
of  which  resemble  the  light  curves  observed  for  W  Ursae 
Majoris  type  stars,  the  orbital  elements  which  they  use  to 
produce  these  similar  light  curves  are  quite  different. 
Unfortunately,  there  is  no  way  to  determine  from  the  two 
synthetic  light  curves  whether  Lucy's  model  and  Wilson's 
model  are  compatible.   The  evidence  of  the  solutions  of 
the  two  synthetic  light  curves  makes  it  clear  that  some 
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caution  is  in  order  in  looking  at  piiblished  solutions  for 
the  light  curves  of  W  Ursae  Majoris  type  stars.   Most  of 
these  solutions  are  based  on  the  geometrical  model  which 
lies  behind  the  Russell-Merrill  method  of  solution  (or, 
equivalently,  other  methods  based  upon  the  same  model) . 
The  solutions  from  this  geometrical  model  (i.e.   the  Russell 
model)  are  not  compatible  with  the  parameters  used  to 
derive  similar  theoretical  light  curves  from  astrophysical 
models  such  as  Lucy's  or  Wilson's.   While  it  is  not  possible 
to  determine  in  this  study  which  of  the  three  models  best 
represent  the  real  situation  in  W  Ursae  Majoris  type  stars, 
there  is  little  doubt  that  there  are  serious  problems  in 
applying  the  Russell  model  to  these  stars. 

One  of  the  problems  encountered  in  the  solutions  of 
both  synthetic  light  curves  was  the  fit  of  the  solutions  at 
the  shoulders  of  the  eclipses.   In  both  cases  the  "obser- 
vations" at  the  shoulders  lie  below  the  theoretical  curves, 
and  the  eclipses  seem  to  be  of  longer  duration  observa- 
tionally  than  predicted  by  the  Russell  model.   According  to 
Bookmyer  "this  particular  discrepancy  between  the  Russell 
model  and  observations  occurs  in  many  W  Ursae  Majoris 
systems"  (16) .   Binnendijk  also  reports  that  "it  has  been 
determined  that  in  many  light  curves  the  observations 
defining  the  shoulders  of  the  eclipse  curves  are  fainter 
than  expected  from  the  Russell  Model"  (1) .   Both  authors 
interpret  this  effect  as  evidence  for  a  permanent  distortion 
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in  the  shapes  of  the  components;  the  facing  hemispheres 
of  the  components  are  more  elongated  than  the  Russell  model 
ellipsoids.   Since  the  synthetic  light  curves  show  this 
feature  which  is  common  to  most  W  Ursae  Majoris  systems, 
perhaps  this  is  further  evidence  that  the  Russell  model 
does  not  apply  to  these  systems. 


CHAPTER  IV 

THE  SYSTEM  OF  44i  BOOTIS 
History 

The  system  of  44i  Bootis  was  discovered  as  a  visual 
binary  by  William  Herschel  in  1781.  On  the  discovery  night 
he  entered  in  his  observing  book: 

Double.   Considerably  unequal.   Both 
white.   With  227  they  seem  almost  to 
touch,  or  at  most  1/4  diameter  of  S. 
asunder;  with  460,  1/2  or  3/4  diam- 
eter of  S.   This  is  a  fine  object  to 
try  a  telescope,  and  a  miniature  of 
Alpha  Geminorum.   Position  29°  54' 
n.   following  (17) . 

The  variation  of  one  or  both  components  was  suspected  by 
several  early  observers  of  the  visual  pair.   These  early 
observations  are  summarized  by  Miss  Agnes  Clerke  as 
follows: 


On  June  16,  1819,  Struve  noted  a 
difference  of  two  magnitudes  between 
the  components;  of  one  invariably 
1822-33,  but  of  only  half  a  magnitude 
1833-38.   Argelander  found  them  ex- 
actly equal,  June  6,  1830.   To  Dawes 
in  April  1841,  the  attendant  star 
seemed  a  shade  brighter  than  its -pri- 
mary, which  was  rated  as  of  fifth 
magnitude.   Duner's  observations  at 
Lund,  1868-75,  confirm  their  relative 
variability,  causing  the  disparity 
between  them  to  range  from  0.4  to 
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1.3  magnitude;  and  he  points  out 
that  they  appeared  to  Herschel  con- 
siderably unequal  in  1781,  but  per- 
fectly matched  in  1787.   Both  stars 
were  yellow  in  1875,  but  the  tint  of 
the  smaller  was  at  times  less  deep 
than  at  others.   Admiral  Smyth 
marked  it  as  "lucid  gray"   in  1842; 
Webb  and  Secchi  respectively  found 
it  blue  in  1850  and  1859;  Webb  and 
Engelmann  reddish  in  1856  and  1865. 
The  principal  star  has  often  been 
considered  as  pure  white.   The 
spectriom  belongs  to  Class  I.   The 
photographic  magnitudes  of  the 
pair,  as  determined  at  Paris  in 
1886,  are  5.3  and  6.   Engelmann 
concluded  the  smaller  component  to 
vary  from  magnitude  5  to  7 ,  the 
larger  from  5  to  6.   They  revolve 
in  a  period  of  261  years,  the  plane 
of  their  orbit  passing  nearly  through 
the  sun.   The  periastron  passage  was 
in  1783.   They  possess  at  least  four 
times  the  solar  luminous  intensity  (18) . 

The  fainter  component  of  the  visual  pair  was  found  to  have 
a  spectrum  resembling  that  of  W  Ursae  Majoris  by  members 
of  the  Mount  Wilson  staff  in  1921.   In  1926  J.  Schilt 
observed  the  star  photographically  and  found  that  this 
component  is  an  eclipsing  variable,  as  suspected,  with  a 
period  a  little  greater  than  a  quarter  of  a  day.   In  1929 
G.  P.  Kuiper  confirmed  the  variability  found  by  Schilt 
from  older  plates  of  the  visual  pair  taken  by  Hertzsprung 
in  1915  and  Munch  in  1922  and  1926.   There  have  been  many 
investigations  of  the  light  variation  of  the  fainter  com- 
ponent of  the  visual  double  since  that  time  (19) . 
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Visual  Binary 

44i  Bootis  is  a  well  known  visual  binary,? 19 09 -ADS 
9494.   The  most  recent  discussion  of  the  orbit  was  by  W.  D. 
Heintz  in  1963.   His  final  elements  as  reported  by  J.  Meeus 
(20)  are 

p   =     246.20  years     a  =     4'.'100 

T   =     2042.00  i   =   +84?50 

e    =       0.360  w   =    237S90 

TT    =       OVOSO  ^      =    228?50 

d 

where  P  is  the  period  of  revolution,  a  is  the  semi -major 
axis  of  the  orbit,  T  is  the  year  of  periastron  passage, 
i  is  the  inclination  of  the  orbit,  e  is  the  eccentricity 
of  the  orbit,  w  is  the  longitude  of  periastron,  tt^  is 
the  dynamical  parallax,  and  ^±s   the  position  angle  of  the 
ascending  node.   Because  of  the  high  inclination  the  ap- 
parent ellipse  is  very  elongated  and  the  angular  distance 
between  the  components  is  less  than  one-half  second  of  arc 
at  this  time  (see  Figure  12) .   The  dynamical  parallax  agrees 
well  with  the  trigonometric  and  spectroscopic  parallaxes 
quoted  by  Binnendijk  (19) . 


TT    =  0'.'076  +  0'.'005 
tr 

7T    =  0'.'073  +  0'.'005 
sp 

We  will  call  the  brighter  component  star  A.   Its  spectral 
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Figure  12.   Visual  Binary  Orbit  of 
44i  Bootis,  (20)  . 
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type  is  F9  (21)  and,  according  to  Binnendijk,  there  are 
indications  that  its  light  shows  a  slight  variability.   The 
fainter  component  of  the  visual  double  star  has  a  spectral 
type  G4  (21)  with  rotationally  broadened  absorption  lines 
like  those  found  in  W  Ursae  Majoris  stars.   We  will  call 
the  fainter  component  (the  eclipsing  variable)  star  B, 
and  its  brighter  and  fainter  components  star  Bl  and  star 
B2  (or  sometimes  star  1  and  star  2)  respectively.   From 
the  visual  orbit  and  a  mean  parallax  of  0'.'076  we  find  for 
the  sum  of  the  masses: 


^A  +  ^Bl  +  ^B2  =  2.59  0  . 


Spectroscopic  Binary 

The  eclipsing  system  is  also  a  spectroscopic  binary. 
D.  M.  Popper  (22)  observed  the  radial  velocity  curves  and 
found  the  following  results : 

Kg^  =  115.4  +1.6  km/sec    K^^   =   231.1  +  1.9  km/sec 

e  =  0  ^2/^Bl  ^  °-^°  -    '^^ 

a^^sin  i  =  0.43  x  10^  km     ^32^^^  i  =  0.85  x  10^  km 
Mg^  sin  i  =  0.77  ©  ^2^^"  ^  ^  0.39© 

Where  K  is  the:  amplitude  of  the  radial  velocity  variation, 
e  is  the  eccentricity  of  the  orbit,  M  is  the  mass  of  the 
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star,  a  is  the  semi -major  axis  of  the  orbit,  and  i  is  the 
inclination  of  the  orbit.   A  more  recent  solution  of  the 
radial  velocity  curves  was  made  by  L.  Binnendijk;  however, 
it  is  based  on  older  observations  made  by  W.  E.  Harper. 
Following  A.  H.  Batten  (23),  Popper's  solution  was  chosen 
because  it  is  based  on  a  greater  number  of  observations. 

From  the  line  intensities  Popper  extimated  that 
the  Iximinosity  ratio  of  the  components  is  0.6  and  most 
certainly  lies  between  the  values  0.5  and  0.85.   Therefore 
in  our  notation  L^  =  0.62  and  L2  =  0.38.   On  the  other 
hand,  Petrie  gives  .06  +  .05  for  the  magnitude  difference 
between  the  components  of  the  spectroscopic  binary  or  0.95 
+0.04  for  the  luminosity  ratio;  thus  L^  =  0.51  and  L2  = 
0.49. 

Eclipsing  yariable: 

A  large  number  of  observations  of  the  light  varia- 
tion of  44i  Bootis  B  have  been  made  since  this  variation 
was  discovered.   These  observations  have  shown  the  presence 
of  strange  deviations  from  the  mean  light  curve.   There 
have  also  been  several  apparently  abrupt  changes  in  the 
period  of  the  system.   There  are  no  published  solutions 
made  from  observations  in  more  than  one  color.   For  these 
reasons,  it  was  decided  that  further  observations  were  in  k 
order. 
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Instriamentation 

The  observations  were  taken  at  the  Rosemary  Hill 
Observatory  of  the  University  of  Florida.   The  telescope 
used  was  a  30-inch  reflector  manufactured  by  the  Tinsley 
Company.   The  telescope  has  a  prime  focal  length  of  120 
inches.   It  has  interchangeable  secondary  mirrors  which 
allow  observing  at  either  the  Newtonian  focus  or  the 
Cassegrain  focus.   A  six-inch  refractor  is  mounted  along- 
side the  reflector  for  use  as  a  finder-scope.   The  photo- 
meter was  mounted  at  the  Cassegrain  focus  with  an  effective 
focal  ratio  of  f/16. 

The  optical  and  mechanical  parts  of  the  photo- 
meter were  built  by  Astromechanics.   As  the  light  rays  pass 
from  the  telescope  into  the  photometer,  they  travel  first 
through  a  small  aperture,  then  through  a  filter,  and 
finally  through  a  Fabry  lens  before  impinging  on  the  photo- 
tube.  The  diaphragm  which  provides  aperture  of  0.33,  0.51, 
0.93,  1.52,  1.98,  2.54,  2.89,  3.90,  4.89,  and  5.85  milli- 
meters, which,  with  a  focal  length  of  480  inches,  corre- 
spond to  angular  aperture  of  5.4,  8.3,  15.2,  25,  32.5, 
41.5,  47,  64,  80,  and  96  seconds  of  arc.   The  32.5-second: 
aperture  was  used  in  these  observations.   A  diagonal  mirror 
may  be  placed  in  the  light  path  following  the  aperture. 
When  the  diagonal  mirror  is  in  the  light  path,  the  light 
rays  are  deflected  to  a  viewing  eyepiece.   The  outline  of 
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the  particular  aperture  in  use  can  be  seen  when  the  eye- 
piece is  in  use,  as  the  aperture  is  illuminated  by  a 
small  bulb.  This  facilitates  centering  the  star  in  the 
aperture.   When  the  diagonal  mirror  is  withdrawn  from  the 
light  path,  the  light  path  is  unobstructed  and  a  signal 
can  be  generated.   When  the  mirror  is  withdrawn  to  take 
a  photoelectric  reading,  the  bulb  is  switched  off  by  a 
micro-switch.   The  starlight  then  passes  through  a  filter. 

The  filters  are  mounted  on  a  wheel  so  that  the 
observer  can  select  one  of  the  ultraviolet  (Corning  9863) , 
blue  (Corning  5030-  and  Schott  GG  13) ,  and  yellow  (Corning 
3389)  filters,  which  correspond  to  the  Johnson-Morgan  UBV 
system  (24) .   The  Fabry  lens  serves  to  keep  the  same  por- 
tion of  the  photocathode  illuminated  despite  the  effect  of 
seeing.   The  filtered  light  rays  impinge  oh  the  cathode  of 
the  EMI  6256B  multiplier  phototube  at  the  end  of  the  light 
path.   A  high  voltage  power  supply  was  used  to  provide  a 
potential  difference  of  1500  volts  between  the  cathode  and 
the  anode  of  the  photomultiplier.   (Unfortunately,  the 
photometer  was  not  wired  properly  by  Astromechanics  and 
there  was  no  potential  applied  to  the  cathode.   Thus  there 
was  no  potential  difference  between  the  cathode  and  the 
first  dynode;  the  whole  1500  volts'  potential  difference  was 
applied  between  the  first  dynode  and  the  anode.   This  fault 
was  not  discovered  until  after  the  observations  had  been 
completed;  however,  the  quality  of  the  observations  does 
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not  seem  to  have  been  affected  by  the  lack  of  potential 
on  the  cathode  of  the  phototube  ..)  The  output  of  the  photo- 
multiplier  tube  is  fed  into  a  DC  amplifier. 

The  DC  amplifier  is  used  to  amplify  the  small 
current  from  the  phototube  to  a  level  detectable  by  the 
recorder.   The  amplifier  was  of  standard  design  with  an 
electric  circuit  similar  to  the  one  described  by  Whitford 
(25) .   The  size  of  the  resistors  is  determined  so  as  to 
change  the  amplification  in  coarse  steps  of  2.5  magnitudes, 
from  0  to  12.5  and  in  fine  steps  of  0.25  magnitudes  from 
0  to  2.50,  where  magnitude  here  refers  to  stellar  magnitude. 
Different  combinations  of  these  resistors  are  used  to  pro- 
duce the  suitable  signal  on  the  recorder  chart.   The  com- 
bination of  resistors  desired  is  selected  with  two  rotary 
switches  mounted  on  the  amplifier.   The  time,  constant  is 
controlled  by  a  third  rotary  switch  which  gives  a  choice  of 
four  time  constants. 

The  output  of  the  DC  amplifier  is  used  to  drive 
the  pen  of  the  dual-channel  recorder  built  by  Honeywell. 
Only  one  channel  of  the  recorder  was  used;  it  has  a  pen 
which  marks  a  continuous  strip-chart  to  correspond  with 
the  signal  current  from  the  amplifier.   The  chart  is  set 
to  run  at  a  rate  of  one  inch  per  minute.   The  observatory 
is  equipped  with  a  high  frequency  receiver  for  receiving 
WWV  time  signals;  the  time  was  recorded  on  the  chart  several 
times  during  the  night.   Readings  of  a  barometer,  a  ther- 
mometer, and  a  hygrometer  were  also  made  several  times  a 
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night  and  recorded  on  the  chart. 

Observations 

The  light  of  the  triple  system  was  measured  by 
centering  the  system  in  a  diaphragm  approximately  thirty 
seconds  of  arc  in  diameter.   Observations  were  obtained 
on  five  nights  in  the  spring  and  summer  of  1969  and  1970. 
Additional  observations  were  made  in  the  spring  of  1971 
for  times  of  minim;am  light.   The  star  HR  5581  was  used  as 
a  comparison  star  throughout.   Its  constancy  had  been 
thoroughly  checked  by  J.  Stebbins  and  C.  M.  Huffer,  and 
later  by  0.  J.  Eggen,  by  comparing  it  with  star  HR  5635 
(26) (see  Table  XXX) .   Binnendijk  (19)  later  reconfirmed 
the  constancy  of  the  comparison  star;  thus,  only  one  check 
star  reading  was  taken  each  night.   The  observations  were 
taken  in  a  regular  sequence:  comparison  star,  sky  de- 
flection, variable  star,  sky,  variable,  comparison,  etc. 
Each  object  was  observed  through  both  blue  and  yellow 
filters, with  individual  readings  lasting  approximately 
thirty  seconds.   Ultraviolet  observations  were  made  on 
only  one  night.   It  was  decided  to  observe  only  in  blue 
and  yellow  light,  since  the  short  period  of  the  system 
requires  observations  closely  spaced  in  time.   The  star 
deflections  were  read  from  the  chart  to  an  accuracy  of 
three  places  (e.g.  54.7)  with  the  sky  plus  dark  current 
deflections  subtracted.   The  time  for  each  observation 
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was  read  from  the  chart  to  the  nearest  5  seconds.   A 
series  of  programs  for  the  IBM  360  computer  was  used  to 
convert  the  raw  data  read  from  the  charts  into  the  final 
light  curve.   The  difference  in  magnitude,  Am,  between  the 
variable  star  and  the  comparison  star  at  the  time  of  a 
variable  star  reading  was  determined  from  the  following 
relation: 

where  i?,(comp.)  is  the  apparent  light  of  the  comparison 
star  and  Jl(var.)  that  of  the  variable  star;  and  this  ratio 
is,  of  course,  equal  to  the  ratio  of  chart-readings  of  the 
two  stars  with  sky  readings  subtracted.   Linear  interpo- 
lation was  used  wherever  values  between  readings  were 
needed.   The  magnitude  differences  together  with  the  times 
of  observation  in  heliocentric  Julian  days  are  given  in 
Table  XXXI.   These  magnitude  differences  were  corrected 
for  first-and  second-order  extinction  and  transformed  to 
the  Johnson-Morgan  UBV  standard  system  in  the  manner 
described  below. 

Reduction  of  Data 

After  observations  have  been  made  of  the  variable 
and  comparison  stars,  it  is  necessary  to  remove  the  effects 
of  the  earth's  atmosphere  from  these  observations.   It  is 
also  most  useful  to  transform  the  data  to  a  standard 
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TABLE  XXXI 
Observations  of  44i  Bootis 


Yellow  Observations 


Hel  JD 

PHASE 

2440000+ 

339.63161 

0.93314 

.63491 

0.94546 

.64.741 

0.99213 

.65048 

0.00359 

.66084 

0.04228 

.66396 

0.05393 

.67444 

0.09306 

.67704 

0.10277 

.68630 

0.13734 

.68879 

0.14664 

.71530 

0.24563 

.71825 

0.25664 

346.61210 

0.99751 

V 6139 5 

0.00442 

.62355 

0.04026 

.62488 

0.04523 

.63183 

0.07118 

.63368 

0.07809 

.64005 

0.10187 

.64173 

0.10814 

.64728 

0.12887 

.64971 

0.13794 

.65579 

0.16064 

.65747 

0.16691 

.66360 

0.18980 

.66545 

0.19671 

.67807 

0.24383 

.67969 

0.24988 

.69861 

0.32053 

.70064 

0.32811 

.70214 

0.33371 

.70804 

0.35574 

.70920 

0.36007 

.71059 

0.36526 

Am 


■0.780 
■0.779 
■0.757 
■0.770 
■0.817 
■0.841; 
■0.878 
■0.881 
•0.870 
■0.875 
-0.862 
■0.907 
■0.765 
■0.767 
■0^802 
■0.814 
-0.826 
-0.845 
-0.857 
-0.865 
■0.882 
-0.882 
-0.885 
■0.893 
-0.901 
-0.901 
-0.899 
-0.909 
-0.878 
-0.881 
-0.887 
-0.878 
-0.881 
-0.867 
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Table  XXXI  continued. 


Hel  JD  PHASE  Am 

2440000+ 


.71806  0.39315 

.71985  0.39983 

.75596  0.53467 

.75868  0.54482 

.77622  0.61031 

.77801  0.61700 

.77980  0.62368 

392.57818  0.62973 

.57997  0.63641 

.59051  0.67577 

.59219  0.68204 

.60029  0.71229 

.60220  0.71942 

.63269  0.83326 

.63455  0.84021 

.64676  0.88580 

.64872  0.89312 

.65781  0.92706 

.65943  0.93311 

.67019  0.97329 

.67210  0.98042 

715.59120  0.99047 

.59462  0.00324 

.60816  0.05380 

.61076  0.06351 

.62262  0.10779 

.62546  0.11840 

.63727  0.16249 

.63975  0.17175 

.65081  0.21305 

.65376  0.22407 

.66481  0.26533 

.66765  0.27593 

.67928  0.31935 

.68188  0.32906 

.69421  0.37510 

.69664  0.38417 

.70735  0.42416 

.71007  0.43432 

.72147  0.47689 

>  .72506  0.49029 

769.59105  0.35379 

; 59325  0.36200 

.59829  0.38082 

.60066  0.38967 


-0. 

863 

-0. 

869 

-0. 

820 

-0. 

826 

-0. 

871 

-0. 

876 

-0. 

881 

-0. 

888 

-0. 

895 

-0. 

935 

-0. 

940 

-0. 

960 

-0. 

970 

-0. 

929 

-0. 

928 

-0. 

903 

-0. 

895 

-0 

834 

-0 

838 

-0 

826 

-0 

832 

-0 

.786 

-0 

.781 

-0 

.859 

-0 

.861 

-0 

.892 

-0 

.885 

-0 

.910 

-0 

.908 

-0 

.928 

-0 

.943 

-0 

.927 

-0 

.926 

-0 

.928 

-0 

.911 

-0 

.883 

-0 

.879 

-0 

.856 

-0 

.855 

-0 

.830 

-0 

.814 

-0 

.944 

-0 

.924 

-0 

.878 

-0 

.871 

Table  XXXI  continued, 
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Hel  JD 

PHASE 

2440000+ 

.60610 

0.40998 

.60847 

0.41883 

.61304 

0.43589 

.61507 

0.44347 

.61906 

0.45837 

.62259 

0.47155 

.62786 

0.49123 

.63035 

0.50053 

.63533 

0.51912 

.63712 

0.52581 

.64175 

0.54309 

.64372 

0.55045 

.64864 

0.56882 

.65089 

0.57722 

.65500 

0.59257 

.65703 

0.60015 

.66142 

0.61654 

.66304 

0.62259 

.66652 

0.63558 

.66860 

0.64335 

.67184 

0.65545 

.67392 

0.66321 

.67769 

0.67729 

.68006 

0.68614 

.68463 

0.70320 

.68642 

0.70989 

.68932 

0.72072 

.69117 

0.72762 

.69539 

0.74338 

.69725 

0.75033 

.70147 

0.76608 

.70326 

0.77277 

.70651 

0.78490 

.70813 

0.79095 

.71154 

0.80368 

.71362 

0.81145 

.71681 

0.82336 

.71872 

0.83049 

.72317 

0.84711 

.72462 

0.85252 

.72867 

0.86764 

.73046 

0.87433 

.73428 

0.88859 

.73590 

0.89464 

.73984 

0.90935 

Am 


-0.842 

-0.859 

-0.835 

-0.828 

-0.801 

-0.826 

-0.811 

-0.796 

-0.815 

-0.821 

-0.827 

-0.828 

-0.864 

-0.879 

-0.883 

-0.890 

-0.892 

-0.916 

-0.890 

-0.907 

-0.930 

-0.917 

-0.937 

-0.931 

-0.942 

-0.967 

-0.931 

-0.924 

-0.894 

-0.908 

-0.907 

-0.916 

-0.913 

-0.926 

-0.928 

-0.907 

-0.906 

-0.921 

-0.917 

-0.913 

-0.919 

-0.914 

-0.870 

-0.844 

-0.860 
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Table  XXXI  continued. 


Hel  JD 

PHASE 

2440000+ 

■•74134 

0.91495 

.74632 

0.93355 

.74782 

0.93915 

.75280 

0.95774 

.75448 

0.96402 

.75969 

0.98347 

.77175 

0.02850 

.78324 

0.07140 

1102.58471 

.58749 

.59392 

.59606 

.60022 

.60283 

.60833 

.61053 

.61463 

.61735 

.62135 

.62366 

.62713 

.62916 

.63333 

.63558 

.64015 

.64230 

.65636 

.65908 

.66568 

.66840 

.67331 

.67557 

.68043 

.68321 

.68715 

.68952 

.69380 

.69571 

.69988 

.70237 

Blue  Observ 

339.62768 

0.91846 

.63590 

0.94915 

Am 


■0.865 
■0.839 
■0.780 
■0.815 
■0.824 
■0.800 
■0.805 
■0.857 
-0.905 
-0.899 
-0.908 
-0.915 
-0.915 
-0.898 
-0.902 
-0.892 
-0.886 
-0.880 
-0.872 
-0.866 
-0.865 
-0.855 
-0.835 
-0.829 
-0.808 
-0.787 
-0.759 
-0.761 
-0.789 
-0.794 
-0.816 
-0.825 
-0.851 
-0.863 
-0.867 
-0.868 
-0.880 
-0.885 
-0.901 
-0.875 


•0.683 
■0.649 


Table  XXXI   continued. 
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Hel  JD 

PHASE 

Am 

2440000+ 

.64672 

0.98956 

-0.647 

.65158 

0.00770 

-0.651 

.65951 

0.03731 

-0.661 

.66466 

0.05654 

-0.688 

.67346 

0.08940 

-0.744 

.eil'il 

0.10624 

-0.734 

.68549 

0.13432 

-0.734 

.68977 

0.15030 

-0.744 

.71449 

0.24260 

-0.748 

.71906 

0.25966 

-0.786 

346.61123 

0.99426 

-0.638 

.61499 

0.00830 

-0.643 

.62292 

0.03791 

-0.671 

.62558 

0.04784 

-0.690 

.63102 

0.06815 

-0.698 

.63443 

0.08089 

-0.720 

.63947 

0.09970 

-0.727 

.64236 

0.11050 

-0.739 

.64664 

0.12648 

-0.750 

.64902 

0.13536 

-0.751 

.65509 

0.15803 

-0.746 

.65816 

0.16949 

-0.755 

.66285 

0.18700 

-0.775 

.66615 

0.19932 

-0.773 

.67732 

0.24103 

-0.762 

.68044 

0.25268 

-0.770 

.71742 

0.39076 

-0.723 

.72072 

0.40308 

-0.720 

.75521 

0.53186 

-0.688 

.75799 

0.54225 

-0.696 

392.57749 

0.62715 

-0.751 

.'58055 

0.63858 

-0.768 

.58987 

0.67338 

-0.800 

.59276 

0.68417 

-0.813 

.59959 

0.70967 

-0.829 

.60283 

0.72177 

-0.829 

.63206 

0.83091 

-0.805 

.63530 

0.84301 

-0.799 

.64606 

0.88319 

-0.768 

.64942 

0.89573 

-0.758 

.65717 

0.92467 

-0.726 

.66007 

0.93550 

-0.699 

.66962 

0.97116 

-0.696 

.67280 

0.98303 

-0.700 

Table  XXXI  continued. 
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Hel  JD 
2440000+ 


PHASE 


Am 


715.59016 
.59566 
.60723 
.61163 
.62181 
.62650 
.63634 
.64051 
.65011 
.65474 
.66400 
.66846 
.67847 
.68275 
.69334 
.69651 
.70659 
.71105 
.72072 
.72592 

769.59152 
.59291 
.59875 
.60020 
.60656 
.60795 
.61357 
.61467 
.61953 
.62080 
.62850 
.62983 
.63573 
.63677 
.64215 
.64331 
.64921 
.65060 
.65541 
.65645 
.66177 
.66258 
.66686 


0.98659 

-0.674 

0.00713 

-0.676 

0.05033 

-0.731 

0.06676 

-0.753 

0.10477 

-0.768 

0.12228 

-0.762 

0.15902 

-0.774 

0.17459 

-0.774 

0.21044 

-0.777 

0.22772 

-0.787 

0.26230 

-0.777 

0.27895 

-0.763 

0.31633 

-0.774 

0.33231 

-0.767 

0.37185 

-0.755 

0.38742 

-0.728 

0.42133 

-0.718 

0.43798 

-0.710 

0.47409 

-0.694 

0.49350 

-0.678 

0.35554 

-0.793 

0.36073 

-0.777 

0.38254 

-0.720 

0.38795 

-0.718 

0.41170 

-0.705 

0.41689 

-0.705 

0.43787 

-0.689 

0.44198 

-0.680 

0.46013 

-0.660 

0.46487 

-0.643 

0.49362 

-0.658 

0.49859 

-0.663 

0.52062 

-0.698 

0.52450 

-0.701 

0.54459 

-0.710 

0.54892 

-0.703 

0.57095 

-0.732 

0.57614 

-0.736 

0.59410 

-0.757 

0.59798 

-0.751 

0.61785 

-0.760 

0.62087 

-0.776 

0.63685 

-0.772 

Table  XXXI  continued. 
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Hel  JD 
2440000+ 


PHASE 


Am,' 


.66820 
.67236 
.67352 
.67821 
.67954 
.68504 
.68614 
.68967 
.69071 
.69580 
.69638 
.70182 
.70292 
.70685 
.70778 
.71200 
.71322 
.71721 
.71831 
.72352 
.72502 
.72902 
.73006 
.73480 
.73561 
.74013 
.74105 
.74673 
.74754 
.75321 
.75408 
.76004 
.77263 
.78284 
1102.58535 
.58686 
.59444 
.59559 
.60080 
.60231 
.60873 
.61006 


0.64186 
0.65739 
0.66172 
0.67923 
0.68420 
0.70473 
0.70884 
0.72202 
0.72591 
0.74491 
0.74708 
0.76739 
0.77150 
0.78617 
0.78964 
0.80540 
0.80996 
0.82485 
0.82898 
0.84841 
0.85402 
0.86895 
0.87283 
0.89053 
0.89356 
0.91043 
0.91387 
0.93508 
0.93810 
0.95927 
0.96252 
0.98478 
0.03179 
0.06991 


-0.777 
-0.802 
-0.772 
-0.779 
-0.797 
-0.813 
-0.819 
-0.799 
-0.797 


-0 


-0 


-0 


-0 


0.777 
0.770 
0.794 
0.781 
0.796 
0.789 
0.788 
0.790 
0.770 
0.762 

791 
0.763 
0.775 

776 
0.723 
0.694 
0.725 
0.715 
0.689 
0.651 
0.717 
0.694 
0.689 
0.692 
0.721 
0.763 
0.751 
0.776 

772 
0.776 
0.770 
0.764 

758 
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Table  XXXI  continued. 


Hel  JD 
2440000+ 


PHASE 


Am 


.61515 
.61660 
.62187 
.62314 
.62766 
.62870 
.63385 
.63506 
.64073 
.64183 
.65694 
.65850 
.66625 
.66770 
.67384 
.67505 
.68101 
.68263 
.68761 
.68900 
.69432 
.69519 
.70051 
.70184 


-0 


-0 


0.739 
0.745 
0.732 
0.733 
0.723 
0.725 
0.695 
0.702 
0.668 
0.654 
0.611 

616 
0.654 
0.665 

681 
0.691 
0.714 
0.716 
0.718 
0.727 
0.737 
0.735 
0.749 

739 


-0 


Ultraviolet  Observations 


715.58808 
.59658 
.60619 
.61250 
.62089 
.62749 
.63547 
.64137 
.64936 
.65590 
.66325 
.66944 
.67760 
.68362 


0.97882 
0.01056 
0.04644 
0.07001 
0.10133 
0.12598 
0.15577 
0.17780 
0.20764 
0.23206 
0.25950 
0.28261 
0.31308 
0.33556 


-0.600 
-0.610 
■0.677 
-0.710 
-0.718 
-0.684 
-0.696 
-0.690 
-0.703 
-0.729 
-0.733 
-0.716 
-0.715 
-0.713 
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Table  XXXI  continued. 


Hel  JD  PHASE  Am 

2440000+ 


.69247  0.36860  -0.828 

.69832  0.39045  -0.689 

.70567  0.41789  -0.648 

.71203  0.44164  -0.644 

.71996  0.47125  -0.642 

.72702  0.49761  -0.634 
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observational  system.   This  transformation  makes  the  results 
independent  of  the  choice  of  telescope-photometer  system, 
thereby  making  direct  comparison  of  data  taken  with  differ- 
ent instrumental  systems  feasible.   The  methods  used  for 
these  reductions  will  closely  correspond  to  those  described 
by  Hardie  (28) . 

Passage  through  the  earth's  atmosphere  has  two 
effects  on  the  light  from  a  star.   First,  it  is  diminished 
and  second,  it  is  reddened.   The  magnitude  of  these  effects 
is  dependent  on  the  amount  of  atmosphere  the  light  passes 
through.   These  extinction  effects  can  be  expressed  as 

m  =  m  -  k'X  -  k"CX  (1) 

o 

where  m  is  the  magnitude  of  the  star  outside  the  atmos- 
phere, m  is  the  magnitude  observed  after  the  light  has 
passed  through  a  path  length  X  of  the  atmosphere,  with  X 
expressed  in  units  of  air  mass  at  the  zenith  of  the  ob- 
server, C  is  the  color  index  of  the  star,  and  k'  and  k" 
we  shall  call  the  first  and  second  order  extinction  coef- 
ficients. 

The  relative  air  mass,  X,  in  units  of  the  air 
mass  at  the  zenith  is  a  function  of  the  secant  of  the 
zenith  distance,  z.   The  value  of  sec  z  may  be  determined 
through  the  relation  . 

sec  z  =  (sincf)  sin6  +  cosp  coscj)  cosh)  (2) 
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in  which  (j)  is  the  observer's  latitude,  while  5  and  h  are 
the  declination  and  hour  angle  of  the  star.   After  sec  z 
is  found  as  above,  then  X  may  be  determined  by  the  formula 
(28) 

X  =  sec  z  -  0.0018167  (sec  z  -  1)  -  0.002875  (sec  z  -1) 
-  0.0008083  (sec  z  -  1) ^  ,  (3) 

To  determine  the  first-order  (or  principal) 
extinction  one  may  measure  a  non-variable  star  over  a 
wide  range  of  air  mass.   Then  the  extinction  coefficient 
may  be  derived  from  a  plot  of  the  measured  magnitude  of 
the  star,  m,  versus  the  air  mass  X.   In  differential  photo- 
metry of  variable  stars,  the  comparison  star  may  serve  as 
the  constant  star  used  to  determine  the  first-order  ex- 
tinction coefficients ,  provided  that  the  observations  cover 
a  sufficiently  large  range  of  air  mass.   Only  on  one  of  the 
five  nights  of  observation  discussed  above  (July  1,  1970) 
were  the  observations  taken  through  a  sufficiently  large 
range  of  air  mass  for  extinction  to  be  determined  accu- 
rately.  The  values  so  determined  (by  a  least  squares  fit 

of  the  data)  were  k'  ,  =  1.15  and  k'   =  1.53.   Plots  of  the 
y  b 

comparison  star  magnitude  versus  air  mass  for  that  night 
are  shown  in  Figure-.13.   For  the  other  four  nights,  it  was 
necessary  to  adopt  the  following  mean  extinction  coeffi- 
cients: k'^  =  0.4,  k'j^  =  0.6,  k'   =  1.0. 

The  second-order  extinction  coefficients  are  a 
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Figure  13.   Comparison  Star  Extinction  for 
July  1,  1970. 


Ill 


measure  of  the  wavelength  dependence  of  the  extinction. 
They  may  be  determined  by  observing  a  pair  of  stars  of  wide- 
ly differing  color  index  over  a  large  range  of  air  mass. 
Then  if  the  two  stars  have  essentially  the  same  position, 
the  differential  measures  of  their  magnitude  and  color  will 
follow  these  relations: 


Am  =  Am  -  k"  ACX  (4) 

o 


AC^  =  AC  -  k"^  ACX  (5) 


Therefore,  plots  of   Am  versus  (ACX)  and  of  AC  versus 
(ACX)  will  give  lines  whose  slopes  are  k"  and  k"  .   These 
second-order  extinction  coefficients  are  substantially 
constant  over  short  periods  of  time  (e.g.  one  observing 
season) .   The  two  stars  of  the  visual  binary  Albireo 
(Beta  Cygni  AB)  were  used  to  determine  the  second-order 
coefficients.   This  system  was  chosen  since  both  stars 
are  relatively  bright,  one  star  is  blue  and  the  other  is 
red,  and  both  stars  are  Johnson-Morgan  standard  stars 
(see  Table  XXXII) .   The  two  stars  were  observed  on  two 
nights  (July  19,  1970  and  August  16,  1970)  through  a 
total  range  of  about  four  air  masses.   When  the  data  was 
reduced  and  plotted  (Figure  14) ,  it  showed  the  presence  of 
a  non-negligible  second-order  extinction  in  the  V  (visible 
or  yellow)  wavelength  band.   This  extinction  is  not  expected; 
however,  it  might  be  attributed  to  the  high  water  vapor 
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TABLE  XXXII 
Standard  Stars 


STAR 


Sp 


o  Tau   G8III 


RA(1970) 


3^23"'lif7 


Dec (1970)   V 


B-V 


U-B 


e  Tau   KOIII  4  26  51.7 

TT"^Ori    F6V  4  48  12.6 

C  Oph    09.5V  16  35  30.2 

a  Oph   A5III  17  33  32.4 

3  Oph   K2III  17  41  59.3 

Y  Lyr   B9III  18  57  49.2 

3  Cyg  A  (KO)  19  29  31.8 

3  Cyg  B  B8V  19  29  31.8 

K  Aql    B0.5III  19  35  16.1 

3  Aql    C8IV  19  53  50.3 


a  Del  B9V 

e  Aqr  AlV 

li  Cap  M2IA 

7 4 Aqr  (B9) 

a  Peg  B9V 

1  Psc  F7V 


20  38  14.6 

20  46  03.1 

21  42  35.2 

22  51  51.8 

23  03  15.8 
23  38  24.3 


8°55.'5  3.59  0.89  0.62 

19  06.9  3.54  1.02  0.88 

6  54.6  3.19  0.45  -0.01 

■10  30.5  2.56  0.02  -0.86 

12  34.8  2.08  0.15  0.10 

4  34.7  2.77  1.16  1.24 
32  38.8  3.25  -0.05  -0.09 
27  53.8  3.07  1.12  0.62 
27  53.8  5.11  -0.10  -0i=32 

■7  05.3  4.96  -0.01  -0.87 

6  19.8  3.71  0.86  0.48 

15  48.3  3.77  -0.06  -0.22 

-9  36.4  3.77  0.01  0.04 

58  38.5  3.99  2.41  2.40 

-11  46.4  5.81  -0.08  -0.32 

15  02.6  2.49  -0.05  -0.06 

5  27.8  4.13  0.51  0.00 
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content  of  the  atmosphere  at  the  time  the  observations 
were  taken  (relative  hiamidity  usually  greater  than  95%)  . 
Since  this  effect  appeared  to  be  an  atmospheric  effect, 
the  k"   term  was  retained  and  the  working  equations  used 
were  modified  to  include  it.   The  second-order  extinction 
coefficients  derived  from  the  data  by  a  least  squaresifit 
were  k"^  =  0.028,  k"^^   =  -0.034,  and  k"^^   =0.022.   Prin- 
cipal extinction  coefficients  were  also  derived  from  the 
observations  of  Albireo  in  order  to  better  determine  mean 
values  for  these  coefficients.   This  data  is  plotted  in 
Figures  15-17. 

Once  the  extinction  coefficients  have  been  deter- 
mined, they  may  be  applied  to  the  observed  magnitudes  in 
the  three  UBV  color  bands  (u,  b,  v)  to  obtain  the  extra- 
atmosphere  values  of  the  stellar  magnitudes  denoted  by 
the  zero  subscripts: 

v^  =  V  -  k'^X  -  k"^(b  -  v)X  (6a) 

b^  =  b  -  k'j^X  -  k"^{h   -  v)X  (6b) 

u  =  u  -  k'^X  -  k"  (u  -  b)X  (6c) 

where  k"^  =  k-^^^  +  k"^  and  k"^  =  k"^^   +  k"^. 

In  order  to  be  able  to  compare  observations  taken 
with  different  equipment,  it  is  necessary  to  transform 
these  quantities  to  a  standard  system  (i.e.  the  Johnson- 
Morgan  UBV  system) .   This  transformation  can  be  represented 
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by  the  following  equations: 

V  =  Vq  +  eCB  -  V)  +  c^  (7a) 

B  -  V  =  y(b  -  v)^  +  ^^  {7b) 

U  -  B  =  ip(u  -  b)^  +  ^j^  (7c) 

where  U,  B,  and  V  are  the  magnitudes  on  the  standard 
system;  e,  u,  and  iJj  are  the  transformation  coefficients; 
and  the  C's  represent  zero-point  differences  between  the 
natural  system  and  the  standard  system. 

One  way  to  determine  these  transformation  coef- 
ficients is  to  observe  a  niunber  of  stars  with  known 
standard  magnitudes  and  color  indices.   These  stars 
should  be  chosen  to  cover  as  wide  a  range  of  color  indices 
as  possible.   The  plots  of  V  -  v^  versus  B  -  V,  B  -  V 
versus  (b  -  v)  ,  and  U  -  B  versus  (u  -  b)^  for  each  star 
will  have  slopes  of  e,   y,  and  ^   respectively.   Seventeen 
stars  chosen  from  the  list  of  Johnson  (29)  were  observed 
and  used  to  determine  the  transformation  coefficients. 
These  stars  with  their  standard  magnitudes  and  colors  are 
listed  in  Table  XXXII.   The  observations  were  corrected 
for  extinction  and  plots  were  made  of  V  -  v  versus  B  -  V, 
(B  -  V)  -  (b  -  v)   versus  B  -  V,  and  (U  -  B)  -  (u  -  b)^ 
versus  U  -  B  (see  Figure  18).   Then  e,  y ,  and  ij;  were 
determined  by  a  least  squares  method.   The  zero-point 
coefficients  were  not  determined.   The  values  of  the  trans- 
formation coefficients  were  e  =  -0.134,  y=  1.072  and  ijj  =  1.''020. 
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Figure  18.   UBV  Transformation  Coefficients. 
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Once  the  extinction  and  trans forroat ion  coefficients 
have  been  determined,  it  is  useful  to  have  a  set  of  working 
equations  which  combine  all  the  reductions.   In  the  study 
of  variable  stars,  it  is  customary  to  select  a  comparison 
star  close  to  the  variable  in  the  sky  and  similar  in  color 
and  magnitude.   By  doing  this  the  differential  measures  are 
less  dependent  on  extinction,  thereby  the  extinction  correc- 
tions are  usually  small.   Considering  equations  (6)  as 
applied  to  both  variable  and  comparison  star  and  expressing 
them  in  differential  magnitudes  and  color  indices: 


Av  =  Av  -  k'   AX  -  k"   A(b  -  v)X        (8a) 
o     ■      V         v 


Ab  =  Ab  -  k',  AX  -  k"   A  (b  -  v)X        (8b) 
o  b        b 


Au  =  Au  -  k'   AX  -  k"   A(u  -  b)X        (8c) 
o  u        u 


where  X  is  the  mean  air  mass  and  third-order  terms  are 
ignored.  Considering  equations  (7)  as  applied  to  both 
stars,  and  rewriting  them  in  terms  of  standard  magnitudes: 

AV  =  Av  +  eijA(b  -  v)i  (9a) 

o  o 

AB  =  Av  +  (y  +  ey)  A  (b  -  v)^  (9b) 

AU   =  Av  +  (u  +  ey)  A(b  -  v)^  +  ipA  (u  -  b)^ 

(9c) 
Equations  (8)  and  (9)  may  be  combined  to  give  the  final 
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working  equations  used  in  the  reductions: 

AV  =  Av  +  (ey(l  -  k"^X)  -  k"^X)A(b  -  v)  -  k'^AX     (10a) 

AB  =  Av  +  (y  +  eyd  -  ^\^^)    "  k"^X)A(b  -  v)  -  k'^^AX 

(10b) 
AU  =  Av  +  (y  +  eyd  -  k"j^^X)  -  k"^X)A(b  -  v) 

+  iP(l  -  k"^j^X)A(u  -  b)  -  k'^AX        (10c) 

The  time  is  the  most  accurately  observed  quantity. 
The  time  in  hours,  minutes,  and  seconds  was  converted  into 
decimals  of  a  day  with  the  day  itself  denoted  by  its  Julian 
Day  number.   These  times  in  geocentric  Julian  day;  must  be 
converted  into  heliocentric  Julian  days  because  of  the 
varying  position  of  the  earth  in  its  orbit  around  the  sun. 
The  light  from  a  star  may  reach  the  earth  as  much  as  eight 
minutes  earlier  and  eight  minutes  later  than  it  reaches  the 
sun,  due  to  the  finite  velocity  of  light.   This  correction 
must  be  made  in  order  not  to  introduce  a  period  change  in 
the  observed  variable  star.   This  correction  is: 

t  =  -O'?005775((cos6  cosa)X  +  (tane  sin6  ; + 
cos6  sina)Y) 

where  a  and  6  are  the  right  ascension  and  declination  of 
the  star,  e  is  the  obliquity  of  the  ecliptic,  X  and  Y  are  the 
rectangular  coordinates  of  the  sun  for  the  given  date,  and 
0cl005775  is  the  time  in  which  light  travels  one  astronomical 
unit  (see  Binnendijk  (30)  for  derivation). 
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Times  of  Minimum  Light  and  the  Period 

The  times  of  minimum  light  were  obtained  from  the 
curves  of  magnitude  difference  as  a  function  of  time  by  the 
tracing  paper  method  (31)  for  the  five  nights  of  observa- 
tion.  An  additional  time  of  minimvim  light  was  obtained 
from  the  observations  made  on  May,  30,  1971  by  the  method 
of  Kwee  and  Van  Woerden  (32).   The  times  of  minimum  light 
are  averages  of  the  values  obtained  with  the  blue  and 
yellow  filters.   They  are  listed  in  Table  XXXIII  along 
with  other  recent  times  of  minimiim  light.   Tabulations  of 
earlier  times  of  minimum  light  are  given  by  Purgathofer 
and  Prochazka  (33)  and  by  Schneller  (34) . 

Variation  in  the  period  caused  by  motion  in  a  visual  binary 
system 

Since  44i  Bootis  B  is  one  component  of  a  visual  bi- 
nary system,  the  observed  times  of  minimum  light  will  vary 
sinusoidally  about  the  true  time  of  minimiam  light  measured 
with  respect  to  the  center  of  gravity  of  the  binary  ".system. 
That  is,  the  observed  time  of  minimum  light  will  differ 
from  the  true  time  with  respect  to  the  center  of  gravity  of 
the  system  by  an  amount  equal  to  the  stellar  light  time 
corresponding  to  the  distance  of  the  variable  from  the 
center  of  gravity  projected  on  the  line  of  sight.   Figure 
19  shows  the  effect  of  this  correction  on  the  selected 
times  of  minimum  light  listed  in  Table  XXXIV.   The  correc- 
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TABLE  XXXIII 
Recent  Times  of  Minimum  Light 


Hel  JD 

E'^ 

0-C" 

2400000+ 

37362.4855 

-4297. 

5 

+0.0009 

7362.6171 

-4297 

-0.0015 

7364.493 

-4290 

-0.0003 

7365.5641 

-4286 

-0.00O4 

7439.6151 

-4009. 

5 

-0.0001 

7443.6322 

-3994. 

5 

-0.0002 

7443.7660 

-3994 

-0.0003 

7449.6577 

-3972 

-0.0005 

7449.7917 

-3971. 

5 

-0.0004 

7472.6889 

-3886 

-0.0013 

7472.8228 

-3885. 

5 

-0.0013 

7740.5042 

-2886 

-0.0003 

7743.4498 

-2875 

-0.0007 

7779.4695 

-2740. 

5 

-0.0020 

7780.4080 

-2737 

-0.0009 

7826.4719 

-2565 

-0.0010 

7835.4445 

-2531. 

5 

-0.0002 

7852.4510 

-2468 

-0.0001 

7853.3867 

-2464. 

5 

-0.0016 

7853.5203 

-2464 

-0.0019 

7855.3959 

-2457 

-0.0010 

7860.4856 

-2438 

+0.0003 

7862.4927 

-2430. 

5 

-0.0012 

7867.4522 

-2412 

+0.0037 

8241.439 

-1015. 

5 

-0.0122 

8540.4649 

101 

-0.0009 

8902.4155 

1452. 

5 

-0.0014 

8910.4510 

1482. 

5 

-0.0003 

8915.4011 

1501 

-0.0048 

8941.3798 

1598 

-0.0040 

9210.3994 

2602. 

5 

-0.0039 

9224.4624 

2655 

-0.0012 

9242.4050 

2722 

-0.0021 

9243.3465 

2725. 

5 

+0.0021 

9245.3520 

2733 

-0.0011 

0-G 
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34 
34 
35 
35 
35 
35 
35 
35 
35 
34 
34 
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34 
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37 
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34 
34 
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Table  XXXIII  continued. 


Hel  JD 

E' 

1 

o-c-^ 

E" 

0-c"  RE 

FEREl 

2400000+ 

9245.4865 

2733. 

5 

-0.0005 

40 

9247.3600 

2740. 

5 

-0.0017 

40 

9247.360 

2740. 

5 

-0.0017 

41 

9248.428 

2744. 

5 

-0.0049 

41 

9248.4310 

2744. 

5 

-0.0019 

40 

9252.4480 

2759. 

5 

-0.0022 

39 

9260.484 

2789. 

5 

-0.0006 

41 

9261.421 

2793 

-0.0009 

41 

9290.3455 

2901 

-0.0004 

39 

9380.328 

3237 

-0.0035 

41 

9391.305 

3278 

-0.0069 

41 

9633.423 

4182 

+0.0070 

41 

9665. 421C 

4301. 

5 

+0.0012 

-2516.5 

+0.0015 

41 

9671. 445c 

4324 

-0.0006 

-2494 

-0.0004 

41 

9673.454c 

4331. 

5 

-0.0002 

-2486.5 

■0.0000 

41 

9760.358C 

4656 

-0.0020 

-2162 

-0.0027 

41 

9922.9250° 

5263 

+0.0017 

-1555 

-0.0008 

42 

9935. 9136^: 

5311. 

5 

+0.0013 

-1506.5 

-0.0013 

42 

9940.8686 

5330 

+0.0018 

-1488 

-0.0009 

43 

9945.4222° 

5347 

+O.0025 

-1471 

-0.0002 

44 

9948.363 

5358 

-0.0026 

-1460 

-0.0054 

45 

9952.9210° 

5375 

+0.0025 

-1443 

-0.0003 

42 

9956.400 

5388 

0.0000 

-1430 

-0.0029 

45 

9959.349 

5399 

+0.0030 

-1419 

+0.0001 

45 

9959.481 

5399. 

5 

+0.0010 

-1418.5 

-0.0018 

45 

9968.453 

5433 

+0.0013 

-1385 

-0.0017 

45 

9977.292 

5466 

+0.0024 

-1352 

-0.0007 

45 

9998.7166 

5546 

+0.0019 

-1272 

-0.0015 

43 

40312.4589 

<6717. 

5 

-0.0003 

-  100.5 

-0.0071 

46 

312.5972 

6718 

+0.0041 

-  100 

-0.0028 

46 

313.4043 

6721 

+0.0078 

-   97 

+0.0009 

46 

316.346 

6732 

+0.0035 

-   86 

-0.0034 

45 

319.296 

6743 

+0.0076 

-   75 

+0.0006 

45 

325.4556° 

6766 

40.0074 

-   52 

+0.0004 

47 

331.346 

6788 

+0.0059 

-  30 

-0.0012 

45 

333.3541° 

6795. 

5 

+0.0054 

-   22.5 

-0.0017 

47 

339.3816° 

6818 

+0.0071 

0 

-0.0001 

47 

339.6493° 

6819 

+0.0070 

1 

-0.0002 

d 

346.6161° 

6845 

+0.0106 

27 

+0.0033 

d 

346.7482° 

6845. 

5 

+0.0088 

27.5 

+0.0015 

d 

363.484 

6908 

+0.0062 

90 

-0.0013 

45 

382.7676 

6980 

+0.0072 

162 

-0.0005 

48 

390.8020 

7010 

+0.0072 

192 

-0.0006 

48 

392.5449° 

7016. 

5 

+0.0093 

198.5 

+0.0015 

d 

Table  XXXIII  continued. 
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Hel  JD 
2400000+ 


0-C 


0-C 


REFERENCES 


392. 

420. 

661. 

675. 

675. 

678. 

679. 

679. 

681. 

681. 

682. 

686. 

697. 

697. 

698. 

698. 

699. 

699. 

700. 

700. 

700. 

700. 

701. 

701. 

714. 

714. 

753. 

769. 

769. 

780. 
1055. 
1055. 
1102. 
1138. 
1139. 
1141. 


6820!? 

7971 

4273 

4875 

6200 

5661 

3684 

5037 

3774 

5119 

4502 

4668 

7184 

8526 

7890 

9213 

7284 

8617 

3943 

6648 

8009 

9324 

7347 

8709 

5891° 
7299  c 

4231 

6286° 

7628° 

471 

3820 

5165 

6556 

406 

345 

4886 


7017 

7122 

8020.5 

8073 

8073.5 

8084.5 

8087.5 

8088 

8095 

8095.5 

8099 

8114 

8156 

8156.5 

8160 

8160.5 

8163.5 

8164 

8166 

8167 

8167.5 

8168 

8171 

8171.5 

8219  ' 

8219.5 

8364 

8424.5 

8425 

8465 

9491.5 

9492 

9668 

9801.5 

9805 

9813 


+0.0125 
+0.0071 
+0.0061 
+0.0071 
+0.0046 
+0.0048 
+0.0036 
+0.0050 
+0.0040 
+0.0046 
+0.0056 
+0.0050 
+0.0084 
+0.0087 
+0.0077 
+0.0061 
+0.0098 
+0.0092 
+0.0061 
+0.0089 
+0.0110 
+0.0086 
+0.0075 
+0.0097 
+0.0068 
+0.0137 
+0.0077 
+0.0104 
+0.0107 
+0.0064 
+0.0060 
+0.0066 
+0.0103 
+0.0075 
+0.0092 
+0.0103 


199 

304 
1202.5 
1255 
1255.5 
1266.5 
1269.5 
1270 
1777 
1277.5 
1281 
1296 
1338 
1338.5 
1342 
1342.5 
1345.5 
1346 
1348 
1349 
1349.5 
1350 
1353 
1353.5 
1401 
1401.5 
1546 
1606.5 
1607 
1647 
2673.5 
2674 
2850 
2983.5 
2987 
2995 


+0.0047 
-0.0011 
-0.0047 
-0.0049 
-0.0063 
-0.0062 
-0.0074 
-0.0060 
-0.0070 
-0.0064 
-0.0055 
-0.0061 
-0.0029 
-0.0026 
-0.0035 
-0.0051 
-0.0015 
-0.0021 
-0.0051 
-0.0025 
-0.0003 
-0.0027 
-0.0038 
-0.0015 
-0.0046 
+0.0023 
-0.0043 
-0.0016 
-0.0013 
-0.0058 
-0.0093 
-0.0087 
-0.0054 
-0.0086 
-0.0070 
-0.0059 


d 
48 
49 
46 
46 
46 
46 
46 
46 
46 
46 
46 
43 
43 
43 
43 
43 
43 
50 
43 
43 
43 
43 
43 
d 
d 
49 
d 
d 
49 
46 
46 
d 
51 
51 
51 


a  Computed  from  min  =  JD  2438513.4166  +  0'?2678143E 

b  Computed  from  min  =  JD  2440339.3817  +  0^26781731E 

c  Used  to  compute  new  light  elements 

d  Present  work 
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Figure  19.   Variation  in  the  Period  of  44i  Bootis  B 
Caused  by  Motion  in  a  Visual  Binary 
System. 
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tions  correspond  to  the  visual  orbital  elements  given 
above  with  the  reasonable  mass  ratio  M  /(M,  +  IVL,  + 
Mg2)  =0.4  with  both  positive  and  negative  inclinations 
of  the  visual  binary  system  being  plotted.   The  computed 
times  of  minimum  light  are  derived  from  an  early  set  of 
light  elements  given  by  Eggen  (26) : 

Min  =  JD2421113.2588  +  0<?26780832  E 

The  plot  shows  several  things.   First  of  all,  the 
sinusoidal  variation  of  the  period  is  overshadowed  by  a 
number  of  apparently  sudden  period  changes;  the  first  one 
occuring  approximately  in  1940.   These  intrinsic  period 
changes  make  it  difficult  to  be  certain  of  the  proper 
choice  for  the  sign  of  the  inclination  of  the  visual  binary 
orbit  and  make  determination  of  the  mass  ratio  of  the  com- 
ponents of  the  visual  binary  unfeasible.   It  appears  likely, 
however,  from  the  observations  before  1940,  that  the  positive 
sign  for  the  inclination  is  the  proper  choice.   Because  the 
light  time  changes  are  nearly  linear  over  the  last  several 
years,  the  uncorrected  times  of  minimum  light  will  be  used 
in  all  further  discussion. 

A  recent  period  change 

The  (0-C)  values  for  all  the  recent  times  of 
minimum  light  (listed  in  Table  XXXIII)  were  calculated 
using  the  most  recent  light  elements  available  as  given 
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by  Pohl  (41) : 

Primary  minimum  =  JD2438573.4166  +  0'?26781430  E 

These  (0-C)  values  are  listed  in  Table  XXXIII  along  with 
the  number  of  elapsed  periods  since  the  given  epoch. 
Figure  20  is  a  plot  of  these  (0-C)  values  as  a  function 
of  the  number  of  elapsed  periods.   From  this  plot  it 
seems  apparent  that  a  sudden  increase  in  the  period  occurred 
sometime  in  1967.   Twenty  times  of  minimum  light  observed 
after  this  time  were  used  to  derive  new  light  elements: 

Primary  minimum  =  JD2440339 . 3817  +  0^26781731  E 

+3  +24 

The  (0-C)  values  and  the  number  of  elapsed  periods 
from  these  light  elements  are  listed  in  Table  XXXIII  and 
plotted  in  Figure  21.   This  new  period  represents  an 
increase  of  about  0.26  seconds  from  Pohl's  period.   The 
phases  of  the  observations  were  now  computed  from  these 
light  elements  and  are  listed  in  Table  XXXI. 

Light  Variations 

When  the  results  of  the  individual  nights  were 
plotted  as  a  function  of  phase  and  combined  (see  Figures 
22  and  23),  several  things  became  apparent.   The  first  was 
the  expected  general  shape  of  the  light  curve;  i.e.  typical, 
of  a  W  Ursae  Majoris  type  variable.   The  second  was  an 
apparent  shift  of  the  zero-point  of  the  magnitudes  from 
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night  to  night;  that  is, the  whole  light  curve  goes  slightly 
up  and  down.   The  third  feature  which  appears  in  the  light 
curve  is  a  slight  depression  shortly  after  maximum  light 
in  the  observations  taken  on  the  night  of  July  1,  1970. 

The  main  feature  in  the  mean  light  curve  is  the 
difference  in  the  height  of  the  two  maxima.   The  secondary 
maximiim  (i.e.  the  one  following  the  shallower  eclipse)  is 
obviously  higher  than  the  primary  maxim\am.   This  has  not 
been  true  in  all  previous  observations  of  the  system.   In 
fact,  the  light  curve  is  subject  to  change  as  is  shown  by 
Figures  24  and  25  which  give  the  mean  light  curves  of  the 
different  authors  who  have  observed  the  system.   The 
earlier  light  curves  up  to  and  including  Binnendijk's  are 
reproduced  from  Binnendijk's  paper  (19).   All  phases  are 
made  zero  at  primary  minimum.   The  zero  points  of  the 
magnitudes  are  adjusted  to  give  the  best  agreement  with 
Binnendijk's  mean  light  curve,  but  the  magnitude  scale 
of  the  authors  was  not  altered. 

Schilt  observed  the  eclipsing  star  with  respect 
to  the  bright  visual  component  photographically.   To 
produce  the  light  curve  shown  in  Figure  24  Binnendijk  made 
normal  points  from  these  long  runs  which  gave  a  good  light 
curve.   Following  Kuiper  he  omitted  three  short  runs  of 
plates  of  poor  quality.   The  light  of  star  1  was  added  in 
the  computation  to  arrive  at  the  combined  light  of  the 
system  using  Schilt 's  value  of  0.76  mag.  for  the  difference 
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Figure  24.   The  Light  Curves  of  Different  Authors, (19) 
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Figure  25.   The  Light  Curves  of  Different  Authors. 
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between  the  visual  components.   The  same  was  done  with  the 
observations  of  the  minima  by  Oosterhoff  (52) .   All  later 
authors  had  to  measure  this  combined  light  because  the 
separation  of  the  stars  became  too  small.   The  light  curve 
by  Stebbins  and  Buffer,  observed  in  1930,  was  computed  by 
Eggen  (26) .   The  normal  points  of  the  light  curve  of  Shapley 
and  Calder  were  only  computed  for  the  observations  taken  in 
1935.   In  Plaut's  light  curve  (53)  Binnendijk  recomputed 
the  normals  around  the  maxima.   Nikonov  gives  light  curves 
in  the  blue  and  infrared;  Binnendijk  reproduces  the  former. 
Normal  points  were  made  from  the  published  graphs  of  Eggen' s 
observations.   This  was  possible  except  for  the  gap  in  the 
figure,  where  the  deviations  from  a  mean  curve  are  as  large 
as  0.04  mag.   Binnendijk  also  omitted  the  deep  and  shifted 
secondary  minimum  which  Eggen  observed  on  March  6,  1947 
from  this  mean  curve.   Binnendijk  computed  the  normal 
points  of  his  own  observations  after  the  zero  points  of 
each  nightVs  observations  had  been  adjusted  to  force  agree- 
ment of  the  maxima.   -Schmidt  and  Schrick  (54)  observed  the 
system  in  three  colors  and  gave  their  observations  in  two 
groups:  Group  I  included  only  good  observations;  and  Group 
II  included  the  remaining  poorer  observations.   The  normal 
points given  here  were  computed  from  the  Group  I  observations 
taken  at  an  effective  wavelength  of  450  my.   Wehlau  and 
Leung  (35)  give  light  curves  in  yellow  and  blue  light;  the 
former  is  reproduced  here.   Catalano  and  Saitta  (36) 
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observed  the  system  in  Ha  light;  their  normal  points 
are  reproduced  here.   Chen  and  Rekenthaler  (55) ,  Brown 
and  Pinnington  (42)  ,  and  the  author  observed  the  system 
in  yellow  and  blue  light.   In  all  cases  normal  points  were 
computed  from  the  yellow  observations. 

An  inspection  of  these  mean  light  curves  shows 
that  both  the  relative  heights  of  the  maxima  and  the  shape 
of  the  light  curve  change.   It  seems  apparent  from  reading 
the  literature  that  the  changes  in  the  shape  of  the  light 
curve  occur  on  a  very  short  time  scale.   Chen  and  Rekenthaler 
(56)  reported  a  variation  in  the  depth  of  the  secondary 
minimum  of  approximately  0.033  mag.  from  one  cycle  to  the 
next,  i.e.  in  a  period  of  less  than  six  and  one-half  hours. 
Other  authors  have  also  reported  night  to  night  changes  in 
the  light  curve.   Eggen,  for  example,  observed  a  deepened 
(0.04  mag)  and  shifted  secondary  minimum  on  the  night  of 
March  16,  1947.   Both  Brown  and  Pinnington  and  the  author 
have  observed  similar  slight  depressions  shortly  after 
secondary  maximum  on  the -nights  of  March  6,  1968  and  July 
1,  1970  respectively  (see  Figure  26) .   These  variations  in 
the  light  curve  and  especially  the  deformity  of  the  light 
curve  at  maximum  light  seem  to  be  present  in  most  of  the 
observations  of  the  system.   Variations  of  light  curves  of 
eclipsing  binaries  of  the  W  Ursae  Majoris  type  are  not  un- 
common.  In  fact,  the  many  irregular  period  changes  of 
the  system  suggest  that  the  system  is  unstable,  and  there- 
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Figure  26.   Deformities  of  the  Light  Curve  of  44i 
Bootis  on  Two  Nights. 
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fore  that  variability  of  the  light  curve  is  not  unexpected. 
In  this  case,  however,  the  picture  is  further  complicated 
by  several  indications  in  earlier  observations  that  the 
brighter  non-eclipsing  component  of  the  visual  pair  is 
slightly  variable.   Untangling  these  two  effects  on  the 
combined  light  will  probably  not  be  possible  until  the 
two  components  of  the  visual  binary  can  be  observed  sepa- 
rately. 

Re  c  t  if  i  cat  ion 


In  order  to  proceed  conveniently  towards  a  solution 
of  the  light  curve,  it  was  necessary  to  convert  the  data 
from  stellar  magnitudes  to  light  values  and  to  free  the 
observations  from  the  light  of  star  A.   There  are  three 
independent  measures  of  the  difference  in  magnitude  Am^g 
between  component  A  and  the  eclipsing  system  B  at  maximum 
light  (21) : 

X  ^^AB  Author 

B  0.76  Schilt 

V  0.63  Kuiper 

V  0.70  Wallenquist 

A  difference  of  0.76  mag.  in  blue  light  and  0.66  mag. 
in  yellow  light  was  used  to  remove  the  light  of  component 
A  from  the  light  curves.   These  values  are  also  consistent 
with  the  spectral  types  of  the  stars  as  reported  by 
Kurpinska  and  Van't  Veer  (21). 
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The  yellow  light  curve  was  chosen  for  the  first 
attempts  at  a  solution.  The  light  in  the  maxima  of  the 
eclipsing  system  was  now  written  in  the  form: 


I  =  A   +  A^cosQ  +  A_cos26  +  A-COS38  +  A.cosf 
O     1  2  J  4 


+  B,sine  +  B2sin2e  +  B^sinSB  +  B^sin4e. 


The  coefficients  were  then  derived  by  using  Merrill's 
graphical  method  as  follows: 

A  =  0.7985 
o 

A^  =-0.0038  B^  =-0.0253 
A2  =-0.0976  B2  =  0.0050 
A^  =  0.0055        ^3  =  ° 

A.  =  0  B.  =  0 

4  4 

The  intensities  were  then  rectified  by  the  formula: 
I  -  A,  cose  -  A-,cos3e  -  B,  sine  -  B^sin2e 

-rll   _   ± f ± ^ 

A   +  A-COs2e 
O     2 

and  the  phase  angles  were  rectified  by  the  formula 

2       2  2 

sin  0  =  sin  e/(l  -  z  cos  e) 

where  z  =  -4Ao/N{A  -  A^)  with  N  =  2.6  corresponding  to  a 
Z  O  z 

limb  darkening  x  =  0.6. 
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Orbital  Elements 

The  rectified  light  curve  shows  two  unequal  minima 

with  light  losses  of  1  -  ^^^  =0.23  and  1  -  ^  =0.13.   After 
^  o  o 

a  coefficient  of  limb  darkening  x  =  0.6  was  assiimed  for  each 
component,  solutions  were  made  with  several  values  of  k,  the 
ratio  of  radii.   As  the  Merrill  nomographs  reveal,  for  shal- 
low minima  of  partial  eclipses  a  large  change  in  k  intro- 
duces only  a  slight  change  in  the  shape  of  the  eclipses. 
This  means  that  for  a  great  part  of  the  partial  eclipse 
portion  of  the  nomographs,  the  shape  (x)  curves  run  nearly 
parallel  to  or  coincident  with  the  depth  line.   If  the 
differences  in  the  depths  of  the  two  eclipses  are  not  great, 
then  the  depth  lines  for  the  two  possible  cases  (the  pri- 
mary eclipse  either  an  occultation  or  a  transit  with  the 
secondary  eclipse  the  other)  will  not  differ  greatly,  be- 
coming coincident  when  the  two  eclipses  are  equally  deep. 
This  indeterminacy,  which  is  illustrated  by  the  nomographs, 
means  that  solutions  must  be  tried  for  a  large  range  of 
values  of  k. 

Solutions  were  made  for  k  =  1.0  and  for  the  two 
cases  of  k  =0.50  (primary  an  occultation  and  primary  a 
transit) .   These  solutions  almost  cover  the  whole  range  of 
k  possible  with  the  assumed  eclipse  depths  and  the  restric- 
tion to  partial  eclipses.   Later  investigation  showed  that 
for  other  possible  values  outside  this  range  (including 
complete  eclipses) ,  the  size  of  one  of  the  components 
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would  exceed  the  size  of  its  Roche  lobe.   Therefore  these 
values  of  k  cover  the  whole  range  of  permitted  solutions. 
Theoretical  light  curves  were  computed  for  the  three 
values  of  k  and  compared  with  the  rectified  observations. 
The  sums  of  the  squares  of  the  mean  deviations  for  all 
three  cases  were  nearly  the  same.   Therefore  the  solution 
for  these  observations  is  indeterminate  over  the  whole 
range  of  values  of  k.   Previous  solutions  in  the  literature 
also  mention  this  problem  of  indeterminacy.   Eggen  found  a 
range  of  indeterminacy  of  k  between  the  two  values  of  0.9, 
although  all  his  solutions  have  one  component  that  exceeds 
or  nearly  exceeds  the  size  of  its  Roche  lobe.   Binnendijk 
found  a  "satisfactory"  solution  with  k  =  0;70  when  he  used 
0.76  for  the  difference  in  magnitude  between  the  eclipsing 
system  and  the  brighter  component  of  the  visual  binary. 
However,  when  he  used  0.63  for  this  difference  in  magnitude, 
he  found  so  large  a  range  of  possible  k's  that  the  solution 
was  indeterminate.   Since  the  light  of  the  brighter  com- 
ponent of  the  visual  binary  cannot  be  measured  more  accu- 
rately at  this  time,  and, in  fact,  this  light  might  even  be 
slightly  variable,  it  does  not  seem  likely  that  further 
observations  at  this  time  will  produce  a  better r determined 
solution. 

There  are  three  solutions  of  photoelectric  light 
curves  of  44i  Bootis  B  in  the  literature.   Eggen' s  solution 
was  made  from  the  unfiltered  photoelectric  light  curve 
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observed  by  Stebbins  and  Huffer  in  1930.   Binnendijk's 
light  curve  was  observed  with  a  yellow  filter  that  does 
not  correspond  to  a  standard  system.   Catalano  and  Saittra 
observed  the  system  in  Ha  light.   These  three  light  curves 
can  be  compared  with  the  author's  yellow  light  curve.   The 
light  curves,  as  illustrated  in  Figure  24  and  Figure  25, 
for  the  four  sets  of  observations  show  somewhat  different 
shapes.   It  is  difficult  to  make  any  direct  comparisons  of 
the  light  curves  since  some  of  the  authors  have  included 
zero-point  adjustments  in  their  data.   It  is  apparent, 
however,  that  differences  in  shape  do  occur  in  the  various 
light  curves.   These  zero-point  adjustments,  and  the 
differing  amounts  of  third  light  make  direct  comparisons 
of  the  various  solutions  somewhat  meaningless.   The  best 
that  can  be  said  is  that  all  the  solutions  suffer  from  the 
lack  of  precise  knowledge  of  the  two  components  of  the 
visual  binary  system  and  from  the  poor  determinacy  of 
solutions  when  the  stars  are  highly  distorted  and  the  geo- 
metrical depth  is  low. 

The  rectification  coefficients  and  the  light  losses 
in  the  rectified  minima  differ  in  the  various  solutions. 
This  is  not  surprising  since  both  the  rectification  coef- 
ficients and  the  light  losses  in  the  rectified  minima  depend 
upon  the  shape  of  the  light  curve  outside  of  eclipse  and 
the  amount  of  third  light  removed  from  the  light  curve.   In 
the  case  of  the  author's  yellow  light  curve,  the  changes  in 
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the  shapes  of  the  light  curve  on  the  five  nights  the 
observations  were  made  add  greatly  to  the  uncertainty  of 
all  the  parameters.   The  light  losses  in  the  rectified 
minima  are  not  well  determined;  however,  this  uncertainty 
in  itself  had  littles  effect  on  the  determinacy  of  the  solu- 
tion.  The  limb  darkening  coefficient  for  both  components 
was  chosen  as  x  =  0.6  following  the  value  determined  for 
the  sun  which  is  of  similar  spectral  type.   AH  three  pre- 
vious solutions  used  a  value  of  limb  darkening  x  =  0.8. 
Changing  the  liirdD  darkening  to  this  value  does  not  signifi- 
cantly affect  the  results  of  the  solution.   There  are  con- 
siderable differences  in  some  of  the  orbital  elements  for 
the  three  solutions  listed  in  Table  XXXV.   While  the  gen- 
eral problem  of  indeterminacy  is  typical  of t shallow  minima 
or  partial  eclipses,  the  very  large  range  of  possible  solu- 
tions may  be  attributed  to  the  large  amount  of  scatter  in 
the  observations.   This  large  amount  of  scatter  makes  it 
virtually  impossible  to  improve  the  determinacy  of  the  solu- 
tion. 

Orbital  elements  were  determined  ^for  these  three 
possible  solutions  of  the  yellow  light  curve.   The  figures 
of  the  components  were  determined  using  the  mass  ratio 
m^/m,  =  0.50  taken  from  Popper's  spectroscopic  elements. 
These  orbital  elements  are  listed  in  Table  XXXV.   The  sizes 
of  the  Roche  lobes  for  the  mass  ratio  ^2/"!^-^   =  0,50  are 
y^  =  0.441,  z^  =  0.414,  ^2   =  0.313,  and  z^   =0.300.   (The 
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TABLE  XXXV 


Solutions  and  Orbital  Elements  with  x  =  0.6,  1  -  It.   =   0.23, 


1  -  i^^^   =  0.13,  0   =  41?0,  and  111^/111^  =0.50 
o  e  z   J. 
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b-  and  c-  axes  of  the  stars  are  in  the  y-  and  z-  directions 
respectively;)   Thus  in  the  solution  where  k  =  0.50  and  the 
primary  eclipse  is  an  occultation,  the  size  of  the  brighter 
star  just  exceeds  the  size  of  its  Roche  lobe.   For  the  other 
case  where  k  =  0.50,  neither  component  exceeds  the  size  of 
its  Roche  lobe.   In  the  case  where  k  =  1.0,   the  fainter 
component  is  larger  in  size  than  its  Roche  lobe.   The  Roche 
lobe  is,  in  effect,  the  limiting  surface  for  a  dynamically 
stable  star  in  a  binary  system  where  the  two  stars  have 
distinct  surfaces.   If  both  components  fill  their  lobes, 
which  may  often  be  the  case  for  W  Ursae  Majoris  systems,  the 
components  will  share  a  common  surface.   The  eclipsing  sys- 
tem of  44i  Bootis  B  is  very  likely  near  this  limit  of  sta- 
bility.  Certainly  the  changes  in  the  period  and  the  light 
curve  indicate  that  the  system  is  not  stable. 

Since  the  solution  of  the  yellow  light  curve  was 
virtually  indeterminate,  no  solution  was  made  from  the 
slightly  poorer  quality  blue  light  curve.   Further  observa- 
tions of  the  system  of  44i  Bootis  for  the  purpose  of  deter- 
mining orbital  elements  will  probably  not  be  productive 
until  the  two  components  of  the  visual  binary  can  be  ob- 
served separately.   When  this  is  possible,  the  problems 
caused  by  the  uncertainty  of  the  difference  of  magnitude  of 
the  two  visual  components  and  the  possible  variability  of 
the  brighter  component  should  be  resolvable.   Then,  it  is 
possible  that  very  accurate  photoelectric  observations 
might  provide  a  better  determined  solution  for  the  orbital 
elements.   Even  under  the  best  circumstances,  however,  the 
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basic  problem  of  indeterminacy  in  shallow  minima  of  partial 
eclipses  will  remain.   Until  the  time  when  the  visual  com- 
ponents can  be  observed  separately,  the  combined  system 
should  be  observed  in  order  to  determine  the  times  of 
minimxom  light,  and  to  detect  from  these  times  of  minimum, 
any  further  period  changes. 


CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 


The  light  curves  of  three  W  Ursae  Ma j oris  sy terns, 
two  generated  from? theoretical  astrophysical  models  and 
one  which  was  observed,  have  been  studied  using  the  Russell- 
Merrill  method  of  solution  of  light  curves  of  eclipsing 
binaries.   The  observed  system,  44i  Bootis  B,  is,  in  many 
ways,  typical  of  W  Ursae  Majoris  systems.   The  light  curve 
of  the  system  changes  with  time  in  an  apparently  irregular 
manner.   The  period  of  the  system  also  shows  irregular 
variations.   These  variations,  which  are  common  in  W  Ursae 
Majoris  systems  may  indicate  that  these  systems  are  not 
dynamically  stable.   In  addition  44i  Bootis  B  is  a  member 
of  a  visual  binary  system. 

Three  major  problems  were  encountered  in  the  solution 
of  the  light  curves.   The  light  curves  of  44i  Bootis  B  and 
the  first  synthetic  system  had  relatively  shallow  minima 
caused  by  partial  eclipses.    This  led  to  a  problem  of 
indeterminacy  in  the  solutions.   The  "observed"  points  in 
the  shoulders  of  both  synthetic  light  curves  fell  below  the 
theoretical  light  curves  predicted  by  the  solutions.   In  all 
three  light  curves  there  was  a  problem  with  third  light.   In 
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the  case  of  44i  Bootis  B,  there  was  a  known  visual  companion 
whose  light  was  included  in  the  observations;  however,  the 
amount  of  this  third  light  was  not  precisely  known.   The 
addition  of  third  light  in  the  solution  of  the  two  synthetic 
light  curves  somewhat  improved  the  fit  of  the  solutions  to 
the  light  curves.   There  was,  however,  no  sound  basis  for 
adding  this  third  light.   In  addition  to  these  problems, 
the  orbital  elements  predicted  by  the  solutions  of  the  two 
synthetic  light  curves  using  the  Russell-Merrill  method  were 
not  at  all  close  to  the  orbital  elements  that  were  used  to 
generate  these  light  curves  from  the  theoretical  astrophysi- 
cal  models.   In  both  cases  the  Russell-Merrill  solutions 
predicted  that  the  stars  were  smaller  and  more  detached  than 
the  actual  stars  used  to  generate  the  light  curves. 

The  results  summarized  above  led  '  to  the  conclusion 
that  the  Russell  model  is  not  compatible  with  the  theoretical 
astrophysical  models  used  to  generate  the  synthetic  light 
curves.   It  is  not  possible  to  determine  from  this  study 
which  of  the  models  more  closely  represents  the  actual  con- 
ditions present  in  W  Ursae  Majoris  systems.   The  solutions 
of  the  light  curves  of  W  Ursae  Majoris  systems  which  are 
in  the  literature  should  be  regarded  as  crude  approximations. 
Better  solutions  may  be  possible  when  the  nature  of  the 
variations  in  the  periods  and  light  curves  of  these  systems 
are  determined.   In  this  light,  it  is  suggested  that  pre- 
cise multicolor  photoelectric  observations  of  a  few  selected 
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systems ,  preferably  those  showing  complete  eclipses ,  carried 
out  as  continuously  as  possible  and  extending  over  a  period 
of  several  years ,  would  be  more  meaningful  than  observations 
of  many  systems  solely  for  the  purpose  of  deriving  orbital 
elements.   These  observations,  together  with  simultaneous 
spectroscopic  observations,  could  add  greatly  to  the  knowl- 
edge of  the  nature  of  the  W  Ursae  Majoris  systems. 
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